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Zur Lösung von Potentialaufgaben mit Hilfe des 
Differenzenverfahrens*) 


Von Johann Schröder in Hannover 


Bei Anwendung des gewöhnlichen Differenzenverfahrens auf Randwertaufgaben der (ebenen oder räumlichen) 
Potentialtheorie erhält man ein lineares Gleichungssystem für die unbekannten ( Näherungs-) Funktionswerte in 
den Rechenpunkten des Gitters. Es wird ein in vielen Fällen anwendbaresVerfahren zur vereinfachten Ermittelung 
dieser Unbekannten angegeben. 


. Applying ihe usual difference method to boundary problems of the (two- or three-dimensional) potential 
theory, the solulion is given by a linear system of equalions satisfied by the (approximate) values in the 
nodal points of the net. The author gives a simpler method of computing the unknowns that is applicable to 
several cases. h 


A l’application de la methode ordinaire de differences a des probl&mes de valeurs ewtrömes de la theorie 
potentielle (bi- et tridimensionnelle) on obtient un sysiöme lin£aire d’Equations pour les valeurs inconnues 
(approximalives) dans les nouds de calcul du reseau. Une methode applicable dans beaucoup de cas & 
l’Evaluation simplifiee de ces valeurs inconnues est indiquee. 


UlpumeueHne OÖ5IYHOTO PA3HOCTHOTO MeTONa K KpaeBbIM 3ANAUYaM (MIOCKOH MIN HPOCTPAHC- 
TBEHHOf) TeOPHH IIOTeHIHAala NMaeT CHTeMy AMHeÄHBX Yp&ABHeHHÄ AA (TPHÖNHKEHHEIX) 
HEH3BCTHEIX 3HAyeHMHH PyHKIIMH B Y3JIOBbIX TOYKAX ceTu. JlaeTca CcIMOCOÖ AA YIPOINeHHOTO 
BbIABJIeHHA 9TUX HEeHSBECTHLIX, IIPHMEHHMBIH BO MHOTHX CIYaAX. 


Für Randwertaufgaben mit Differentialgleichungen der Gestalt AU—cU+p=0 
(e = const.) läßt sich die Rechenarbeit bei Benutzung des gewöhnlichen Differenzenverfahrens 
und Verwendung eines Quadrat- oder Sechseck- bzw. Würfelgitters in vielen Fällen wesentlich 
herabdrücken. Bei geeigneter Numerierung der Rechenpunkte des verwendeten Gitters und 
der zugehörigen (Näherungs-) Funktionswerte bekommt man ein Gleichungssystem einer 
speziellen Form. Gleichungssysteme solcher Gestalt werden allgemein im Abschnitt 1 behan- 
delt. Die dort erhaltenen Ergebnisse werden im Abschnitt 2 auf Randwertaufgaben verschie- 
dener Art mit einer Differentialgleichung der oben beschriebenen Form angewendet. Der 
Abschnitt 3 bringt Beispiele. 


1. Vereinfachte Lösung von linearen Gleiehungssystemen einer speziellen Form 
A. Aufgabenstellung 


Gegeben sei ein lineares Gleichungssystem 
n 
— ou + I bu tr = 0 Ve Den ta) 
Dt 


für die Unbekannten w;. Darin seien entweder &, b;, und r; gegebene Konstanten (Fall I), oder 
es liege eine Eigenwertaufgabe vor (Fall II), bei der & gleich einem Parameter x ist, die b;x 
bekannte Polynome dieses Parameters sind und alle Zahlen r; verschwinden. Vorausgesetzt 
wird, daß eseine natürliche Zahl 1 < n gibt, mit der 


een lundmteLl,k>l . ..n ..%, 0) 


B. Beduktion der Anzahlder Unbekannten 


a) Ist B die n-reihige quadratische Matrix mit den Elementen b;; — B = (bi,) — und 
sind u und r die n-dimensionalen Vektoren mit den Komponenten u; bzw. n — u= (u), 
r = (r;) —, so lautet (1) in Matrizenschreibweise 


ei Se ee rl Ba Er c)? 
Mit E, (» positiv, ganz) bezeichnen wir die v-reihige Einheitsmatrix. 


*) Über diese Arbeit trug der Verfasser auf der GaMM-Tagung 1953in Aachen vor. 
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Wir fassen die ersten l und die restlichen n — 1 Komponenten eines n-dimensionalen 
Vektors zu je einem neuen |- bzw. (n — l)-dimensionalen Vektor zusammen, z. B. 


“\ u r [ "+1 
Ug u+2 en 15 IR 17+2 

Ur | WU g ’ er== ZE® = a ’ 
U Un r, In 


v s SR 
und schreiben dann in leicht verständlicher Symbolik u = (2) ‚= (4) und entsprechend bei 


anderen Vektoren. Analog dazu denken wir uns auch n-reihige quadratische Matrizen in vier 
Teilmatrizen zerlegt. Von den beiden aus den ersten I bzw. letzten n — 1 Zeilen gebildeten 
Matrizen werden jeweils die ersten / und letzten n — I Spalten zu einer neuen Matrix zusam- 
mengefaßt, z.B. 


b1,1+1 b1,142° "' bı,n bi+1,1bı+1,2°*  biyn.ı 
Br b2,1+41 ba,1ı+2 ° * ' be,n | Bi bi+2,1 bi+2,2 ° * * bir, ı 
b,,1+1 Duıro bin bn,1 bn,2 u ba, 
Symbolisch werde geschrieben B = 1 2 und entsprechend bei anderen Matrizen. 
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Mit dieser Schreibweise geht (3) über in 
lie ER ENAN 


ee ee 
Br-en tree ET: 


E,| B 
b) Multipliziert man die Gl. (4) !) mit der Matrix € = ü er 
n 


lee 


(- ME + BB)v+o=0 mt DeasH Bi Tre Eee 
BR u Hal EEE 


was gleichbedeutend ist mit 


) ‚ so ergibt sich 
0 


- B, 


..(6) 
Gh 


Jede Lösung u = der Gl. (4)genügt auch den Beziehungen (7). Im Falle det C = 


—a+0,d.h.x +0, erfüllt auch umgekehrt jede Lösung von (7) die Gl. (4). It «=0, 
so enthält (7) keine Bedingung für w. 

Das Gleichungssystem (7a)!) besteht aus / Gleichungen für die Unbekannten 4,, %, . . . ‚ur. 
Da bei diesem System also die Anzahl der Unbekannten (wie der Gleichungen) gegenüber dem 
Ausgangssystem (3) von n auf I reduziert ist, nennen wir es das ‚‚reduzierte System‘“ Hat man 
im Falle +0 die Unbekannten «,, 4... , u,aus (7a) berechnet, so ergeben sich die restlichen 
Werte 4741, %4+2, .. ., 4, sofort einzeln aus (7b). Die von (4) ausgehende Ermittelung des Systems 
(7a) erfordert u.a. eine Matrizenmultiplikation, die Bildung der Matrix B, B,; bei vielen der 
im Abschnitt 2 besprochenen Anwendungen läßt sich das System jedoch direkt aufstellen. 


c) Ist (1) eine Eigenwertaufgabe (Fall Il, «= x, b; =b; (2), M— 0), so haben die 
Gleichungen (3) und (7a) die spezielle Gestalt 
(ED) l ER Eee 
bzw. 

(— mE, + BB)u—= 0 (9), 


1) Beziehungen dieser Art bezeichnen wir jenach Bedarf als Gleichungen — nämlich Matrizengleichungen — 
oder Gleichungssysteme — nämlich Systeme gewöhnlicher Gleichungen. 


; Fi. ri 
3 
fi 


u ae ER 


& w u 


folgenden Beziehung 


i hle welche die a. 8. einen Rt aen Rektor besitzt, sind die 
Be "Nullstellen der Funktion g (x) = det (— x E„+ B). Auf Grund der au 


7) = = det Sir, 


= x En-—ı 7 det ( — E07, | Nn— | 
= det 1 »®E,+ BıB,) : (— „ni 


En-ı B, 


ii ) - det ei 


det (— “BE, +B) = (— Dar ®!det (— a2; +B,B,) . . (10) 


erhält man diese Zahlen auch mit Hilfe der Determinante des reduzierten en . Für 
eine Zahl x +0 gibt es wegen (5) (dort speziell & = x, s = t = 0 gesetzt) keinen von Null ver- 


v 
schiedenen Lösungsvektor u = bei dem einer der Vektoren v oder w verschwindet. 


Genügt jedoch u = (2) für«—=0den Gin. (5), so werden diese auch von den Vektoren (5) 
und (2) erfüllt. 


Hängen die Zahlen b;, nicht von x ab, so sind es die ‚‚charakteristischen Zahlen“ x, der 
Matrix B, für welche die Gl. (8) einen von Null verschiedenen Lösungsvektor — dann ‚‚Eigen- 
vektor‘ der Matrix B genannt — besitzt. Aus (10) folgt, daß dann mit x auch — x charakte- 


ristische Zahl ist. Und wegen (5) ist S Eigenvektor zur Zahl — x, falls (2 ) solcher zur 
Zahl x ist. ” 


C. Anwendung von Iterationsverfahren 


a) Es liege der Fall I (&, b;2,r; bekannte Zahlen) vor, und es gelte «+0 und 
det (—a&E„-+ B) +0. Die Systeme (3) und (7) besitzen dann die gleiche eindeutige Lösung 


u= |,,); man kann versuchen, sie iterativ zu ermitteln. Beim Ausgehen von (7) wird man ein 


Iterationsverfahren natürlich nur aufdaserste Teilsystem (7a)für die Unbekannten u,, %..., % 


anwenden und die restlichen Unbekannten anschließend aus (7b) berechnen. Bei dem im 


Folgenden durchgeführten Vergleich von Iterationsverfahren für die Gleichungssysteme (3) 
und (7a) werden wir nur das Verhalten der Näherungen für die ersten ! Komponenten 
des Lösungsvektors untersuchen und dieses als charakteristisch für die Schnelligkeit der 
Konvergenz ansehen. Die verschiedene Zahl der Gleichungen bei (3) und (7a) bleibt unberück- 
sichtigt. ; 


k 
Für die zwei Folgen der n-dimensionalen Vektoren u* = E und der !-dimensionalen 


0 
Vektoren (k=0,1, 2,...), welche von einem Vektor u® -(&) bzw. dem Vektor g’ = v 
ausgehend nach der Vorschrift % 


’ RED Ne A er er (11) 
ZW. 
i BE 2A ae a ee ee ee (12) 
gebildet werden, gilt die Beziehung 
De ee U (13). 


(11) ist die Vorschrift für das Iterationsverfahren in Gesamtschritten zum Gleichungssystem 
(3). Herrscht für dieses Verfahren Konvergenz, so konvergiert das durch (12) charakterisierte 
Iterationsverfahren also doppelt so schnell in dem durch die Gl. (13) beschriebenen Sinne. 
Man darf annehmen, daß das Iterationsverfahren in Gesamtschritten für das reduzierte 
System gewöhnlich noch etwas schneller konvergiert, da die Matrix B, B, im Normalfall 
von Null verschiedene Glieder in der Hauptdiagonalen besitzt, die in (12) als Faktoren 
der Komponenten g* auftreten, während sie beim Gesamtschrittverfahren mit zu den 


Faktoren der g#t! geschlagen werden. 


Meistens verwendet man nicht das Iterationsverfahren in Gesamtschritten, sondern das 
in Einzelschritten. Auch bei diesem ist jedoch im allgemeinen eine erhebliche Konvergenz- 


k 
verbesserung zu erwarten, Sind nämlich «* = (2) (k=0,1,2,...) die mit Hilfe der Vor- 
schrift 


ziS 


u 


einen Eigenvektor, wenn zwei verschiedene charakteristische Zahlen der Matrix auch ver- 


Me EURER 


es 2 ii, 2 pi 
Schröder, Lösung von 


Eu 


gebildeten Vektoren — diese Vorschrift entspricht de 
(3) bzw. (5) bei Beginn mit der (1 + 1)-ten Zeile —, so gilty = 
Einzelschritten für das Ausgangssystem konvergiert also (bei B te 
so schnell wie das durch (12) beschriebene Verfahren, d.h.im: Igemeinen langsan 
Verfahren in Gesamtschritten für das reduzierte System und damit gewöhnlich erheblich 
samer als das Einzelschrittverfahren für dies System ?). er: f 
b) Es werde jetzt der Fall II (Eigenwertaufgabe) unter der zusätzlichen j 
betrachtet, daß die Zahlen b;, und damit auch die Matrix B von x nicht abhängen un diersTe 
charakteristischen Zahlen dieser Matrix sämtlich reell ausfallen. gr E DE R. 
Zur Berechnung eines Eigenvektors und einer zugehörigen charakteristischen Zahl einer . 
Matrix benutzt man häufig ein Iterationsverfahren, bei welchem man eine Folge von Vektoren 
durch wiederholte Multiplikation eines geeignet gewählten Ausgangsvektors mit der Matrix 
bildet®). Diese Vektoren konvergieren, mit geeigneten Faktoren multipliziert, immer dann gegen 


schiedene Beträge haben ®). Für die Matrix B ist jedoch diese Voraussetzung nicht erfüllt (es 
seien denn alle charakteristischen Zahlen gleich Null), da mit x auch gleichzeitig —x charakte- 
ristische Zahl von Bist 5). In diesem Falle zeigen die genannten Vektoren, welche also hier von 


: | EN 
einem Vektor u — (3) ausgehend nach der Vorschrift +! = Buf(k=0,1,2,...)gebldet 


werden, gewöhnlich ein alternierendes Verhalten. Die charakteristischen Zahlen der im redu- 
zierten System auftretenden Matrix B, B, sind aber gerade Quadrate von charakteristischen 
Zahlen der Matrix B, also auf Grund unserer Voraussetzung sämtlich positiv, so daß für die mit 
einem Vektor g® nach der Vorschrift +!= B,B,g* (k=0,1,2,...) gebildeten Vektoren 
immer Konvergenz in dem oben angegebenen Sinne herrscht. Ist g° = v", so gilt hier ent- 
sprechend (13) q* = v?*. In diesem Falle erhält man also bei Iteration mit der Matrix des 
Ausgangssystems (8) nach 2% Schritten die gleichen Näherungen für die Unbekannten «,, 
ugs... , u Wie bei Iteration mit der Matrix des reduzierten Systems nach k Schritten. 


2. Anwendung auf Randwertaufgaben 


A. Aufgabenstellung und Bezeichnungen 


a) Es sei eine Randwertaufgabe folgender Art gegeben. @seiein beschränktes Gebiet in 
einer %, y-Ebene [einem x, y, z-Raum] mit einem Rande I’. Gesucht ist eine bestimmte Stetig- 
keitsforderungen erfüllende Funktion U (z, y)[U (x, y,2)], welche in @ der Differential- 
gleichung 

AU—cU+p=0 


mit AU=U.+ Up [Ur + Uyy + U] genügt und auf T’ gewisse gegebene lineare Rand- 
bedingungen erfüllt, über welche unten noch etwas genauere Aussagen gemacht werden. Dabei 
sei entweder peine gegebene Funktion der Veränderlichen z, y [®, y, 2] und c eine gegebene 
reelle Konstante (Fall I) oder p=0 und c darstellbar als c= c«,— A mit einer gegebenen 
reellen Konstanten c, und einem Eigenwert A (Fall II, Eigenwertaufgabe). 

Bei Anwendung des gewöhnlichen Differenzenverfahrens auf dieses Problem ergibt sich 
im folgenden ein System von linearen Gleichungen (18) und (20). Wir wollen annehmen, daß 
sowohl die gegebene Aufgabe als auch dieses Gleichungssystem lösbar sind (im Falle I ein- 
deutig) und daß das gewöhnliche Differenzenverfahren auf Näherungslösungen für die gestellte 
Aufgabe führt‘). Zur vereinfachten Ermittelung dieser Näherungslösungen, d. h. der Lösungen 
des genannten Gleichungssystems, kann man das im folgenden beschriebene Verfahren be- 
nutzen. Dies Verfahren wird dargestellt an ebenen Problemen bei Benutzung eines quadra- 
tischen Gitters. Die folgenden Ausführungen gelten jedoch auch für Sechseckgitter und ebenso 
für räumliche Probleme bei Benutzung eines Würfelgitters. Man hat nur jeweils die entspre- 
chenden, in eckigen Klammern angegebenen Abweichungen zu benutzen und ausdrücklich nur 
für quadratisches Gitter gemachte Aussagen wegzulassen. 

b) Näherungswerte für U und A werden mit kleinen Buchstaben u bzw. A bezeichnet. 
Zur Lösung der gegebenen Aufgabe mit Hilfe des gewöhnlichen Differenzenverfahrens werde 

°) Vgl. [1], J. Schröder, Eine Bemerkung z. Konvergenz d. Iterationsverf -f. lin. Gleichgssyst., Arch. 
Math. IV (1953), 8. 322—326. 


®) Vgl. [2), L.Colla tz, Eigenwertaufgaben mit tecbnischenAnwendungen, Leipzig1949, S.294 ff. u.S.309f. 
4) Siehe [2], S. 309 f. 
5) Allerdings kann man durch eine Verschiebung des Spektrums Konvergenz erreichen, indem man also 
statt B eine Matrix B — u E„ mit einer geeigneten reellen Zabl u betrachtet. 
°) Zu Theorie und Anwendung des Differenzenverfahrens siehe etwa [2], S. 336 ff. oder [3], L. Collatz 
Numerische Behandlung von Differentialgleichungen. Berlin 1951, S. 279. ; 
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7 8 " Summe der (Näherungs-) Funktions 
y; werte in den Nachbarpunkten 

mit 

s—=ß+ch? und ß=4 [ß = 3beim Sechseck-, ß — 6 beim Würfelgitter] . . . (14b). 


Diese Punkte sollen ‚‚innere Gitterpunkte‘ genannt werden. Es gebe deren mindestens zwei. 


Wir setzen voraus, daß die inneren Gitterpunkte „zusammenhängend“ sind; d.h. zu je zwei 


nichtbenachbarten inneren Gitterpunkten Pund P*solles innere Gitterpunkte P’, P”,..., Pw 


(u >1) geben derart, daß in der Folge P, P', P",..., Pw, P* je zwei nebeneinanderstehend 


genannte Punkte Nachbarpunkte sind. Alle Nachbarpunkte von inneren Gitterpunkten, welche 


nicht selbst innere Gitterpunkte sind, heißen ‚‚Randpunkte“. In jedem Randpunkt werden die. 


gegebenen Randbedingungen herangezogen, um eine weitere Gleichung für die Funktionswerte 
zu erhalten. Es werde vorausgesetzt, daß jede dieser Gleichungen die Gestalt 


Li . . . = . 1 
1 | inearkombination aus den (Näherungs-) Funktionswerten in +y=0..(15) 


den benachbarten inneren Gitterpunkten und Randpunkten 


hat, wo uder Wert in dem betreffenden Randpunkt und yeine bekannte Konstante (im Falle II 
gleich Null) ist. Diein der Linearkombination auftretenden Koeffizienten können im Falle II 
Polynome in Asein. Ineinem Randpunkt, für welchen die Linearkombination in (15) verschwin- 
det, ist der Wert vbekannt, und dieser kann in die Gln.(14) und die übrigen Gln. (15) eingesetzt 
werden. Die inneren Gitterpunkte und die Randpunkte, in denen w nicht von vornherein 
bekannt ist, diejenigen Gitterpunkte also, in denen die Funktionswerte berechnet werden sollen, 
nennen wir „‚Rechenpunkte‘‘ ?). Ihre Zahl sein. Wir bezeichnen sie mit ; = 1,2,..., n), 
u; seien die zugehörigen Näherungsfunktionswerte. Allgemein soll ein Funktionswert in einem 
Gitterpunkt wie dieser Punkt gekennzeichnet werden; win P’, u*in P* usw. 

c) Für ebene Probleme kommt man b»i Benutzung eines quadratischen Gitters z. B. in recht allgemeinen 
Fällen auf Gleichungen der Form (15) 8). 


In Bild 1sei Pein Randpunkt, P’und P’ bedeuten benachbarte Rechenpunkte. In dem Schnittpunkt P* 
der (in einer Umgebung von P* als glatt vorausgesetzten) Randkurve /’mit einer durch den Punkt P gehenden 
Geraden sei die Randbedingung . 


AUGSBURG N 1. Saas era (16) 


gegeben, worin U*% diein P* genommene Ableitung der Funktion Uin Richtung der inneren Normalen von I’be- 
deutet. Durch einen Index u wird im folgenden die Ableitung in Richtung PP’’’ gekennzeichnet. Im Falle A = 0 
sei @ so gewählt, daß diese Richtung mit der Richtung der inneren Normalen in P* zusammenfällt. Für A= 0 
braucht dies.nicht der Fall zu sein, im allgemeinen wird man dann = 0 verwenden. Es soll eine lineare Beziehung 
zwischen «, w’ und w’’ aufgestellt werden. Beilinearer Interpolation erhält man 
sing uw’ + cos ou” 
sing + cos 
Setzt man andererseits als Näherung 
w’— U*+b ur und u= U*—a Urs 


ZZ 
U 


so ist 
A— Bb AB. 
re 
Im Falle A — Bb = Oerhält man durch Elimination von w’”’ 
A+Basinpw +cospu’ td _ 
A—Bb sinp+cosp A—Bb 
welches eine Gleichung der Form (15) ist. Liegt P zwischen P’’ und P*, so ist a negativ zu nehmen, liegt P’’’ 


zwischen P* und P, hat man b negativ anzusetzen. Auch der Winkel p kann negativ sein oder den Winkel /2 
überschreiten. 


w’ +0= AU +BUr+C=0. 


ak ae, Bild 1 


?) Wir benutzen bier die bei L.Collatz[3]S. 281 verwendeten Bezeichnungen in geringfügiger Abwand- 
lung. ‘Dort werden zur Darstellung der Theorie des Differenzenverfahrens bei der ersten Randwertaufgabe alle 
Randpunkte als Rechenpunkte bezeichnet. 

8) Zur Berücksichtigung der Randbedingungen siehe auch etwa [4], F.B. Hildebrand, Methods of 
applied Mathematics, New York 1952, S. 301 ff. 


HRp=0%. 2.2... (ide) 


a) Die Punkte P; werden so 
zwei benachbarten Punkte zur gleichen S 


kurz „„A-Punkte“ — und Pıyı, Pı4z,..., Pn die Punkte 
sind (im allgemeinen wählt man die Bezeichnung der Sorten so, daß I = 


b) Die Differenzengleichung (14a) in einem inneren Gitterpunkt P; 
läßt sich schreiben als 


ee 


‚der: 


E 


n. ' iR En SUN, B 
— u + IS buumtn=0.... 0.200000. li) S 

k=1 k S. k Ay > ER 

mit : ! { . | AZ ri K: 2% i 

Are 1 falls P,dem Punkt P; benachbartist  ,, (19a 3 Be: 

Be 0 sonst ng 
E- % | nei + h? pi; a ee le ee Se a ET De er ee (19b), a : & 
er t; = Summe der Funktionswerte in den Nachbarpunkten Ri: 
= von P;, die keine Rechenpunkte sind, = 
& > pr = pP (m, y) Ip (@ Yu %)]- WERE “ 
E Jede Gl. (15) in einem Randpunkt P; hat die Gestalt - 5 


= + 2 pr pi — 0 N N j ande (20), 


‘ [3 
wobei die @;, höchstens für solche Werte k von Null verschieden sind, für welche P; ein P; 
benachbarter Rechenpunkt ist. Der Rechenpunkt P;, in welchem die Gl. (18) bzw. (20) auf- 
gestellt ist, soll ,,Aufpunkt‘ dieser Gleichung genannt werden. rs 
c) Die Gln. (18) in den inneren Gitterpunkten und die mit & multiplizierten Gln. (20) 
bilden ein System, das sich in der Gestalt (1) 


n 5 - - 
— aut Ib nt tr; — 0 Me er = 8) 
k=1 E 


schreiben läßt (für einen Randpunkt P; ist b;,; = % 9;, und r; = ay;) und für welches infolge 
der Art der Numerierung der Punkte P; auch die Voraussetzung (2) erfüllt ist, so daß wir 
die Ergebnisse aus Abschnitt 1 benutzen können (der Buchstabe « hat dort eine andere 
Bedeutung!). Außer «; kommen ja in den Gln. (18) und (20) jeweils nur Funktionswerte in 
Punkten vor, welche nicht zur gleichen Sorte wie P; gehören. Im Falle x + O ist dieses System 
äquivalent dem System der Gln. (18) und (20). In diesem Fall können daher nach Abschnitt 1 
auch das reduzierte System (7a) und das bezüglich der Unbekannten +1, Wa, - - . , %, ent- 
koppelte System (7b) zur Berechnung der Werte u; dienen. Diese lauten in Summenschreib- 


weise 
1 
ort I Bir unt 0=0 a PA EN eg (22a), 

bzw. % 

— 0 +2 br rn 0 = li +1,12... SS . (22b), 
wenn ;, die Elemente der Matrix B,B, und o; die Komponenten des Vektors o sind, 

n n 
ir => b;;b;x und eg =ıar-+ Sn b;;r; Of OS eg Sl fen (23). 
j=1+1 j=l+H 


Zur Ermittelung der Werte u;hätte manim Falle x + Oalso die Größen ß;, und g;aus (23) 
zu bestimmen (d.h. im wesentlichen die Matrix B, und den Vektor t mit der Matrix B, zu 
multiplizieren), dann das reduzierte System (22a) für die Funktionswerte in den A-Punkten zu 
lösen und anschließend die restlichen Unbekannten aus (22b) zu berechnen. Die GlIn. (22a) 
lassen sich in vielen Fällen jedoch auch direkt, ohne Benutzung von (23) aufstellen. Darauf 
wird unter Ü — E eingegangen. 

Zur Kontrolle für die im System (22) auftretenden Koeffizienten kann man die Zahlen 


I I 
= 2 Pr -entaei=l2...,)) und b=Sbutn G=1+1,1+2...,n) (24) 


ER Au a nun im TE WITTEN Le Ausai ı EE: 2% 
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berechnen, mit ihnen muß gelten 


da wegen (23) : 
l h n tz n 

L = 3 Pr -ın to = D Bu 2 but n;) = 2 b,b;. 

k=1 = j-IH1 - 


je \= 
‚d) Liegt eine Eigenwertaufgabe vor (Fall II, c= «, — A), so ist & ein zu bestimmender 
Parameter, welcher dann «heiße. Erist mit dem (Näherungs-) Eigenwert Adurch die Gleichung 
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verknüpft. Gesucht sind Zahlen x, für welche die (n-reihige) Determinante D des Systems 
der Gln. (18) und (20) verschwindet. Da die GIn. (20) mit « multipliziert wurden, um sie auf 
die Gestalt (21) zu bringen, ist diese Determinante gleich ®) «" det (— x6;, + b; x) , wenn m 
die Anzahl der zuden Punkten P; gehörigen Randpunkte ist, und damit nach Abschnitt 1 


gleich der mit (— 1)*-!,%-2!== multiplizierten (l-reihigen) Determinante des reduzierten 
Systems (22a), 


D=(— 177 „nam det (— 226;5+ Pix) EEE cn al a (27). / 


e,) Vielfach kann man auf Grund von Symmetrieeigenschaften des Gleichungssystems 
(14) (15) bzw. der gestellten Aufgabe aussagen, daß die Werte u und u* in gewissen Rechen- 
punkten P und P* (im Falle II wenigstens für bestimmte Lösungen) gleich sind, u = u*. 
Gehören diese Punkte zur gleichen Sorte, so lassen sich derartige Beziehungen zur Herab- 
setzung der Zahl der Gleichungen und Unbekannten in (22a) verwenden 1°). 

Wir denken uns die Rechenpunkte wie bisher in zwei Sorten Aund B geteilt, bezeichnen 
deren Zahl aber jetzt mit N. Auch die Numerierung geschehe in etwas anderer Weise. 
Unter diesen Rechenpunkten gebe es n Punkte — P,, P,,..., P„ genannt — mit der Eigen- 
schaft, daß für jeden Rechenpunkt P mit dem Funktionswert u eine natürliche Zahl? <n 
existiert, derart daß (im Falle II jedenfalls für die Lösung, welche man berechnen will) 


ist und P und P; zur gleichen Sorte gehören. Die in irgendeinem Rechenpunkt aufgestellte 
Differenzengleichung (14) bzw. (15) stellt dann eine Beziehung zwischen den Werten uy 
(k=1,2,...,n)dar. Es wird ferner vorausgesetzt, daß diese Beziehungen für die Punkte P 
und P; identisch sind. Gewöhnlich wird man eben auf Grund der Identität dieser Gleichungen 
auf (28)schließen, bzw. auf Grund von Symmetrieeigenschaften der gestellten Aufgabe, welche 
die Identität dieser Gleichungen zur Folge hat. Die Punkte P; werden nun unter sich so nu- 
meriert, daß mit geeigneter Zahl! P,P,,..., P, A-Punkte und die übrigen Punkte B-Punkte 
sind. 

e,) Man hat dann nur in den jetzt mit P; bezeichneten Punkten die Gln. (14) bzw. (15) 
unter Benutzung der Beziehungen (28) aufzustellen. Diese Gleichungen haben wieder die 
Gestalt (18) bzw. (20), wobei die Konstanten b;, und @;, auch hier nur für solche Werte k von 
Null verschieden sind, für welche P; und P; zu verschiedenen Sorten gehören. Es gilt daher 
auch hier das oben unter c) und d) Gesagte. 


Dieser ste Randwertaufgabe.bei passendem Grieer 


a) Es sei U auf I’ gegeben. Ferner werde angenommen, daß alle Randpunkte auf 7’ 
liegen, daß das Gitter auf den Bereich ‚‚passe‘“. Dann sind alle Rechenpunkte innere Gitter- 
punkte. In diesem Falle lassen sich die Zahlen ß;, und g; in (22a) leicht direkt ermitteln 
(vgl. Bild 3, wir verwenden nun wieder die alten Bezeichnungen): 


ß;x = Anzahl derjenigen Rechenpunkte P;, die sowohl Nachbarpunkte von P;als 
auch von P; sind (kurz: Anzahl der gemeinsamen Nachbarn); für «= k 
soll dies bedeuten, daß ß;; die Anzahl der P; benachbarten Rechenpunkte ist . . (29a), 


0; = ar; + Ir, wobei über alle jsummiert wird, für die P;und P;benachbartsind . . (29b). 
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10) Bei quadratischem Gitter kann man auf die im folgenden beschriebene Weise Symmetriegeraden berück- 
sichtigen, welche mit Seiten oder Diagonalen von Quadraten des Gitters zusammenfallen. Symmetriegeraden, 
welche die Seiten von Gitterquadraten jeweils in der Mitte zwischen zwei benachbarten Gitterpunkten senkrecht 
schneiden, sind in anderer Weise zu berücksichtigen. Man kann sie als Rand auffassen, auf dem die Ableitung UR 
nach derinneren Normalen verschwindet, und die Ergebnisse aus Abschnitt D benutzen. Geraden, bezüglich derer 
die Lösung antimetrisch ist, lassen sich als Rand behandeln, auf welchem U= (0. 
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b) Die Näherungsgleichung (14) in 


renzenstern veranschaulichen, bei quadra 
jede der GIn. (22a) kann man einen zugehörigen > 
hängt, welche Nachbarpunkte von P; Randpunkte 
in der Tabelle angegebenen Sterne 1 bis 5 vor. 
c) Da unter den Rechenpunkten kein Ran ist, hängen die 
Falle II (Eigenwertaufgabe) nicht von x ab. Die gesuchten Zahlen x sin 
ristische Zahlen der symmetrischen Matrix B = (b;x) reell. Als Wurzeln < Gl 
(D wie in (27) mit m = 0) liegen sie daher symmetrisch zu x = 0 auf der reel 


r 


erhält man aus einem zur charakteristischen Zahl x gehörigen Lösungssystem (%:) ein zur Zahl 
— xgehöriges, indem man die Werte in den B-Punkten mit —1 multipliziert. Ist x=0 charak- 
& teristische Zahl, so bekommt man die zugehörigen Lösungen gewöhnlich nicht alle mit Hilfe der 
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d) ©) us 
ae m dom Kreis Bie- Bias. Aufgabe: Gesucht ist eine Funktion U (z, 9), welche in dem bei a akizxierten Geblet der Poiential- 


itter- gleichung A U = 0 genügt und auf dem Rande die angegebenen Werte 1 bzw. — 1 annimmt. Die Rechen- 
a Fe an Diese punkte ?,, P,,..., P, deseingezeichneten Gitters sind durch ihre Nummern 1,2,...,?7 gekennzeichnet. 
h? pdar Die Zahlen ß,, B;n 0; dj: dj, 7; sind bei c in der bei d erklärten Anordnung aufgeführt 


Gin. (22). Im allgemeinen geben derartige Lösungen und ebenso die zugehörige Zahl A aber 
auch nur schlechte Näherungen für U bzw. Aab. Zu x= 0 gehörige Lösungssysteme (%;) 
setzen sich ja nach Abschnitt 1 linear aus solchen zusammen, bei denen die Funktionswerte 
entweder in den A-Punkten oder in den B-Punkten gleich Null sind. Die Form solcher Lösun- 
gen hängt also sehr eng mit der Gestalt des gewählten Gitters zusammen. Man kann derartige 
Lösungen aber oft auch leicht direkt angeben unter Benutzung der Tatsache, daß bei einer 
Lösung, für die alle Werte in den Punkten einer Sorte verschwinden, in jedem Punkt dieser 
Sorte die Summe der Funktionswerte in den benachbarten Punkten gleich Null sein muß. 

d) Im Falle I konvergiert für das Gleichungssystem (21), falls e > 0, sowohl das Itera- 
tionsverfahren in Gesamtschriften als auch das in Einzelschritten 4). Nach Abschnitt 1, C 
kommt man bei Iteration des reduzierten Systems aber im allgemeinen mit erheblich weniger 
Iterationsschritten zur gleichen Genauigkeit für die Funktionswerte in den A-Punkten 2): 
Im Falle II konvergiert das übliche Iterationsverfahren 2) bei Anwendung auf das reduzierte 
System (22a). Hier erreicht man nach %k Schritten die gleichen Näherungen für die Werte in 
den A-Punkten wie bei Anwendung auf das System (21) mit 2 % Schritten. 


' 41) Siehe [3], 288 ff. und [5], L. Collatz, Über die Konvergenzkriterien bei Iterationsverfahren für 
lineare Gleichungssysteme, Math. Z. 53 (1950), S. 149— 161. 

12) Bei Lösung des reduzierten Systems hat man also neben dem Vorteil der geringeren Zabl der Unbe- 
kannten den der schnelleren Konvergenz der Iterationsverfahren. Allerdings muß man dazu sagen, daß die Matrix 
des reduzierten Systems gewöhnlich stärker besetzt ist, d. h. sie enthält pro Zeile mehr von Null verschiedene Rle- 
mente. 
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des passenden Gitters teile man dien Punkte P; in zwei Sorten 4A und Bund 

jere sie wie unter B, a) beschrieben. 

n B-Punkten Pi Pırz, ..., Pu stelle man Haie "gewöhnlichen Diener 

22 gleichungen (22b) auf (unter Benutzung von ceb) und (19) oder des zugehörigen Diffe- 

-  - renzensternes). 

»; Man ermittele die Zahlen r; mit (19b) auch für i=1,2,...,lund stelle die Gln. (22) 
auf (mit Hilfe von (14b) und (29) oder zugehöriger Differenzensterne, wie sie für 
den Fall quadratischen Gitters in der Tabelle angegeben sind. 


c) Zur Kontrolle berechne man mit (24) die Zahlen ß; und b; und stelle fest, ob damit (1) = 


gilt. Außerdem kann man (30) nachprüfen.) 
3) Man berechne: 
a) nur bei Eigenwertaufgaben («= „): die Eigenwerte A aus den Gleichungen D = O und 
(26), wo D’ durch (27) mit m = 0 gegeben ist, 
== und füra #0: 
b) die Unbekannten «,, us, .... u durch Lösung des reduzierten Systems (22a), 
c) dierestlichen Unbekannten aus den Gin. (22b) oder als 1/& - Summe der Funktions werte 
in den Nachbarpunkten. 


f) Besitzt die Aufgabe Symmetrieeigenschaften von der unter B, e,) beschriebenen Art 14), 
so wollen wir unter P;@=1,2,..., n) jetzt die dort so bezeichneten Punkte verstehen. Das 
reduzierte System (22a) erhält man dann, was hier nicht bewiesen werden soll, indem man in 
den Punkten’B.; PP.’ .P; Differenzengleichungen mit Hilfe der in der Tabelle angegebenen 


Sterne aufstellt [bei Verwendung eines Sechseck- oder Würfelgitters mit Hilfe entsprechender, 


durch Überlagerung gewöhnlicher Differenzensterne leicht zu SU ENTED Sterne] und dabei 
die Symmetriebeziehungen der Art (28) benutzt. 

Die in (22a) auftretenden Konstanten ß;, und o; lassen sich aber auch wieder durch 
(29) entsprechende Gleichungen ermitteln. Die Punkte P; bezeichnen wir jetzt auch als 
P,@=132,...,n). Jeder nicht zu diesen Punkten gehörige Rechenpunkt P bekommt als 
ersten Index eine Zahl :, mit der die Gln. (28) und die dort weiter genannten Voraussetzungen 
gelten, und wird von dem Punkt P; und den evtl. noch vorhandenen Punkten mit dem gleichen 
ersten Index durch einen anderen zweiten Index unterschieden. Dann gilt 


- Pr = = Pi, UNO O0 KERLE. 2.0 (32), 


wobei man die Zahlen ß;,,, und g,, aus den Gin. (29) erhält, wenn man dort 2, 7 und k durch % 
7, bzw. k, ersetzt und r;, = r; benutzt. 

Die oben unter c) und d) gemachten Aussagen gelten auch hier bei folgender Abänderung 
unter c). Die Matrix B = (b;,) braucht jetzt nicht symmetrisch zu sein. Ihre charakteristischen 
Zahlen sind jedoch auch hier reell, da das Gleichungssystem (21) hier aus einem Gleichungs- 
system mit symmetrischer Matrix unter Benutzung von Gleichungen der Art (28) entsteht. 
Auch die Beschreibung des Rechenganges unter e) ist hier zu verwenden, wenn man die in 
Klammern stehenden Angaben wegläßt. Es gilt jetzt nicht mehr allgemein (30). Jedoch ist 
Pi, von Null verschieden, falls fz;; esist. Anstelle von (31) kann man zur Kontrolle die folgende 
Gleichung benutzen: 


ß: = & b;, mit b,;, — b,, wobei über alle 5, summiert wird, 
» "für die P,, und P; benachbart sind. 


D. Krummlinige Ränder, andere Randbedingungen 


a) Es werde jetzt zugelassen, daß sich unter den Rechenpunkten auch Randpunkte 
befinden. Wirsetzen jedoch voraus, daß jeder Randpunkt P; nureinem inneren Gitterpunkt 


_ P,benachbart ist und daß die gegebenen Randbedingungen in P; auf eine Gl. (20) der speziellen 


Gestalt 


führen, worin u; der Funktionswert in eben dem P,; benachbarten inneren Gitterpunkt P; ist. 
Dieser innere .Gitterpunkt heiße „dem Randpunkt P; zugeordnet“. Die Gl. (17) nimmt für 
9 = Oeine derartige Gestalt an. Der Winkel 9 = 0 ist dort zu verwenden, falls in (16) A=0 
ist, d.h. die erste Randwertaufgabe vorliegt, oder falls die Randkurve I’ die Seite PP’ senk- 


recht schneidet. 


14) Vgl. auch 10). 


„ 


- 


250 Schröder, Lösung von Potentialaufgaben mit Hilfe des Differenzenverfahrens rn 


Auch in diesem Falle lassen sich die Gln. (22a) durch Differenzensterne veranschaulichen, 
bei quadratischem Gitter kann man die Sterne in der Tabelle verwenden. Ist für jeden Rand- 
punkt P; die in der zugehörigen Gleichung auftretende Konstante g; gleich 1, so kann man die 
Zahlen ß;, und oe; auch wieder verhältnismäßig einfach durch (29) entsprechende Gleichungen 
erhalten (eine Gl. der Form (33) mit 9; = 1 entsteht aus (17) mit B=0undp = 0). Auf deren 
Angabe sei hier jedoch verzichtet. Der Rechengang verläuft ähnlich wie unterC, e) beschrieben. 


Differenzensterntabelle 


Der Stern ist (mit —= 4 + ch?) zu lesen als 


Mw + 2[us +up + we + ua] + [we +ur + ug + ur] 


1 +ah2p; +h2[pa + PB + pc + po]l= 0 
=, —02 +4 
Mw +2[up + wel + [ua +ua tur t ug + ur] 
2 toayatok?p +h?[pg + pc + ppl= 0 
M=—02+3-+09y4 
Mw+2uH+lu+tustUu +urltoalya + vB] 
3 +oh?p +h?[pco + pp]l= 0 
M= —-0o?+2+talgaı + ps] 
Mu;+ [ua +up tue +ua +ur + ur] 
4 +alya tyol+tah:pi + Rh? [ps + pp]= 0 
= —0?+2 +alpa+ pc] 
Mutl[luwe +ua +url +alyat+tys+vye] 
bi} +ah?p, +h?2pp= 0 
M=—a+1l+oa[lyı +tgs + go] 
Mu; t+agilue +ua + ur] +02 y; 
6 +9%hpp= 0 
a 
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3. Zwei Beispiele 


Aufgabe 1. Gesucht ist eine Funktion U(z, y), die im Innern des skizzierten, von zwei 
konzentrischen Quadraten begrenzten Bereiches der Potentialgleichung AU = 0 genügt, auf 
deminneren Rande den Wert U = 1 undauf dem äußeren Rande den Wert U = O annimmt °°). 
Diese Funktion ist offenbar symmetrisch zu jeder der gestrichelt gezeichneten Linien. Wir 
werden hier zwei Quadratgitter verwenden, die ‚„‚passen‘ und symmetrisch zu diesen (gestri- 
E chelten) Symmetriegeraden liegen und bei denen außerdem zwei zu einer solchen Geraden 
symmetrische Rechenpunkte zur gleichen Sorte gehören. Dann brauchen wir nur die Werte u 
in denjenigen Rechenpunkten zu berechnen, welche im Innern des dicht schraffierten Bereiches 
bzw. aufseinem Rande liegen und können nach B, eund (, fdiese mit P;= P;(e=12,...,n) 
bezeichnen und das reduzierte System mit Hilfe der Sterne in der Tabelle oder mit Hilfe der 
GlIn. (32) aufstellen. 
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Bild 6 
2 

Als erstes verwenden wir ein Gitter mit der Maschenweite (Seitenlänge) Ah a a 

(s. Bild 6a) und drei A- und fünf B-Punkten Pf =1,2,...,nn=8) Esista=4. Die 


übrigen im Gleichungssystem (22) auftretenden Zahlen sind in Spalte 1 der Zahlentafel 1 auf- 
15) Dies ist eine Übungsaufgabe aus [3], S. 353, 356. 


ee (ooh) a 
punkten zu berechnen. 


- Genauere Werte wird man nn B 

schenweite (Seitenlänge) h = YAaı 

n = 15) erhalten. Es ist wieder «—=4. Die Za 
om u reduzierten Systems wird Jetzt mit Fulfe 


el gan Gitters ae dan gleichen Nummern ae kann man als i 


Spalte 2 angegebenen benutzen. Dreimalige Iteration liefert Näherunge 


Differenzen ur — us (Spalte 4b) sind für fast alle « dem Betrage nach etwa 
Differenzen u? — ul (Spalte da) und von gleichem Vorzeichen 16)“ Man wi 
sich die folgenden Differenzen ut! — u: angenähert wie ‚Glieder der Form con: 


Zahlentafell e > b wu 
1 En 3a : Fr 
Pil- Br Inia| . 8 
Fe ; SEN uU ik 
b;| bir ri N i 
12) 5-22. |°0.'3 0,3966 411102222 1 
121° 472.020. 4 0,4655 24000440 1 
9) 404 011 0,2155 20424002 0 
1) 001. |0\ 0,0539 00222000 0 
21.1.0110 0,1530 2:0,2.1%3.12.0%1 0 
2| 1100 0,2155 22001410 0 
4, 2.1.0.2 0,7155 22000141 6 
TU 0,5022 40202024 2 = 
4a 4b 4c 4d 5a 5b 
Bear — 105 2 5 & 
wu A; Re = 
607 313 0,39217 | 0,39256 | 0,1530 | 0,1530 
439 281 0,45109 | 0,45170 | 0,1539 | 0,1534 
478 239 0,20154 | 0,20170 | 0,15057 | 0,15213 
Ba 82 | 0,04982 | 0,04993 | 0,14912 | 0,15109 
333 145 0,14767 | 0,14782 | 0,14842 | 0,15023 
295 143- 0,21089 | 0,21117 | 0,14810 | 0,14959 
289 144 0,69743 | 0,69765 | 0,14795 | 0,14912 
402 193 0,47182 | 0,47205 | 0,14788 | 0,14878 


Zahlentafel2 


1 2a 2b 2c 2d 2e 3 4 
u; ob. Vorz.| ug; UT) ul; i 
IR u Us; ti Sr 
« ug zunt.Vorz.| Ur; Us; “ > uı a 
8,632 1 0 1 0 0 1 8,592 
16 0,8328 1 —1,133 0 0 0,9247 18,746 
26,520 0,8828 —l —1,133 0 0 0,9247 30,56 
32 +0,3022 | 0,5 | 1,654 1 1 0,3162 
32 +0,6446 0,5| 0,19 | —ı1 | —1 0,6547 
37,480 +0,6446 | +0,5| 0,194 1-1 0,6510 
48 0,6446 +0,5| F0,194 | —1 1 0,6966 
95,368 0,9869 *0,5| +1,266 l lt 1,0696 


geometrischen Reihe verhalten und entsprechend extrapolieren, d.h. die letzte Differenz 
w — ui noch einmal zu u} addieren, und bekommt auf diese Weise Werte 7 %; (Spalte Ac), welche 
gegenüber den durch weitere Iteration und Extrapolation sehr schnell zu erhaltenden Werten «; 


16) Nur für i—= 4ist eine starke Abweichung davon festzustellen. Die Betrachtung der Zahlen ur, u und u? 
legt jedoch nahe, für = 4 nicht anders zu verfahren, als für die anderen Werte i. n ne 


Je 
n,so bekommt man z. B. für us ee dis 
ungen w #=0,1,...,7). Bei Iteration. des Ausgangssystems (2 
"bei Beginn mit der %ten M + 1)-ten) Zeile (vgl. 1, C, a) erhält man unter 
g der gleichen Ausgangswerte uU (= 1,2,. 
gleiche die Zahlen i in den Spalten 5a und b miteinander und mit der Zahl u.. 


Au f 8 abe 2. (Eigenwertaufgabe). Das Gebiet @ ist das gleiche wie in Aufgabe il 
Gesucht sind zu den gestrichelt gezeichneten Geraden symmetrische, auf dem gesamten Rande 
verschwindende Funktionen U (a, y) und Zahlen A, derart daß AU+-AU=0in G. Wir 
verwenden dieselben Gitter wie in Aufgabe 1 und erhalten dann auch die.gleichen Zah- 
len Br und b;,. Jtztista ma =4— Ali, n = —=0. 


Beim ersten Gitter mit der Maschenweite B= x a erhält man als Nullstellen = sich 
aus (27)ergebenden Determinante D= — x2(2?— 4)(x?— 9x2-+4) die Zahlen x, — V: + 5.165, 


= 2, = 3518, N KB=—h Muh, B——H ee auf 


2 
die in Spalte 1 der ne 2 stehenden Eigenwerte A, führen. Zugehörige Lösungen u,; sind 
in den Spalten 2a — e aufgeführt (das obere Vorzeichen gehört dort jeweils zu den Werten «,) 
mit kleinerem »). Für v #4 oder 5 erhält man diese Lösungen mit Hilfe des Systems (22. 


Fürv=4 und v=5 ergeben sie sich leicht unter Benutzung der Tatsache, daß die Werte in 


den A-Punkten selbst und die Summe der Nachbarwerte in den A-Punkten gleich Null sein 
muß. 


Für das zweite Gitter mit der Maschenweite h — ni bekommt man von den Werten 1; in 


Spalte 2a ausgehend durch Iteration ®) R reduzierten Systems nacheinander Zahlen «1,, u1: 
u; und us; (Spalte 3). Die Quotienten — — „liegen alle zwischen 11,981 und 12,002, innerhalb 


welcher SeHtranken nach dem nen ns (die größte charakteristische Zahl der 
Matrix B, B,) x ?]iegen muß. Mit dem Näherungswert x? = 11,99 wird Ay a? = 8,592 (Spalte 4). 
Die Zahlen «,; in Spalte 2b sind Lösungen des reduzierten Systems für 4 — 8, woraus /, a = 
18,746 folgt (Spalte 4). Exakt ist A,a=2n?= 19,739. Die Zahl x3 wurde mit Hilfe des 
Iterationsverfahrens in sechs Schritten bei Ausgehen von den Werten us;in Spalte 2c zu 4,36 
mit einer Genauigkeit von etwa 2% berechnet (A, a®in Spalte 4). Dabei wurden vor dem ersten 
und sechsten Iterationsschritt die Komponenten des 0-ten bzw. 5-ten Näherungsvektors in 
Richtung der in den Spalten 3 und 2b angegebenen Vektoren (uj;) und (u,;) abgezogen 19). 
Der Wert für x? und die Genauigkeitsaussage ergaben sich mit Hilfe der Schwarzschen 
Konstanten und Quotienten sowie einer bei L. Collatz angegebenen Fehler-Abschätzung 
unter Benutzung der Annahme, daß A, > 32 19) 20). 


17) Diese Zahlen sind auch bei [3], S. 356 angegeben. 
18) Siehe [2],8.291 oder [6],L.Collat z, Einschließungssatz für die charakteristischen Zahlen von.Matrizen. 
Math. Z. 48 (1942), S. 221 — 226. 


19) Als „inneres Produkt‘‘ zweier !- eenaler VekioIeR (u;) und (v;) wurde dabei nicht das gewöhnliche 
[; 


innere Produkt Pr u;v;, sondern der Ausdruck y- & e;u;v; benutzt, wo «; die Anzahl der Rechenpunkte 
= 1 
P;, (i= 1,2, -- +) bedeutet (= 5, = = 4 = 8, ,= = 4= = 4) und y die positive Konstante 1/4 ist. 


20) Siehe [2], S. 298 ff. 


Eingegangen am 14. 6. 1953. 


‚ 8) die Zahlen % in Spalte 5b. Man nver- 
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Graphische Lösung der Hauptaufgabe der 


. . E . F \ ur: e 5 r #)1) ” = 
E  Luftphotogrammetrie im Sinne der Ausgleichsrechnung”? hi 
x . N Se * 2 . = 
wer Von Josef Krames in Wien ra en ° 
= Für die Wiederherstellung der gegenseitigen Lage zweier aus verschiedenen Lufistandpunkten auf- 
pr en Geländebilder ed Er ach graphisches Verfahren abgeleitet. Dieses gründet sich auf 


Pr die geometrische Theorie jener „‚gefährlichen‘‘ Geländeflächen, für welche diese Aufgabe in der Praxis nur 
ungenau bestimmt ist. ’ u Er 


= A handsome graphic method is deduced for the restitution of the relative position of two aerial pholo- 
. graphs taken from two different centres. This method bases upon the geometrical theory of the so-called „eriti- 
cal“ ground-surfaces for which the considered problem practically allows only an inewact solution. 


Une simple methode graphique est &tablie pour la restitution de la position relative de deux vues aeriennes 

5, du terrain, prises de deux centres differents. Cette methode est deduite de la_theorie geomötrique des surfaces 

Be topographiques particulieres, dites ‚‚critiques“, pour lesquelles, en pratique, le probleme envisage ne permei 
qu’une solution inezacte. - 


? B HacToaımei CTATbe BEIBONUTCA IPocToA rpadmyeckuä CIIOCOO AA BOCCTAHOBJIEHHA B3AUM- 
HOTO PACHOJIOKEHHA ABYX a8PocheMOR MeCTHOCTH, 3ACHATEIX C ABYX PasIMYHHX IIOSUNHÄ. ITOT 
CIOCOD OCHOBEIBAETCH HA TEOMETPHYeCKOH TeOPHH TeX „KPHTHYECKHX“ ILIOIMANEH MeCTHOCTH, IA 
KOTOPHX 9T& 3a1a4ya HA IIPAKTHKE UMeeT JIHIIb HETOYHOE OIPpeeleHHe. 


Die Hauptaufgabe der Luftphotogrammetrie besteht bekanntlich darin, je zwei nach- 
einander aus der Luft aufgenommene Geländebilder in jene gegenseitige Lage zu bringen, die 
sie bei den Aufnahmen inne hatten. Der Abstand der beiden Aufnamezentren, die Basis- 
strecke, ist dabei auf eine geeignete Länge 5b zu verkleinern. Die Herstellung dieser „gegen- 
seitigen Orientierung‘‘ zweier Luftbilder oder der mit ihnen verbundenen Zielstrahlbündel 
bildet in erster Linie ein geometrisches Problem, und für die praktischen Anwendungen über- 
dies eine Aufgabe der Ausgleichsrechnung. 

Wir setzen zunächst voraus, daß die beiden Zielstrahlbündel im Orientierungsgerät 
mathematisch exakt wiederhergestellt seien. Ihre gegenseitige Orientierung ist (abgesehen 
von gewissen ‚„‚gefährlichen‘‘ Geländeformen s. unten) er- 
reicht, sobald im Raume jeder Strahl des einen Bündels von 
seinem zugeordneten im anderen exakt geschnitten wird. 
Das aus den Schnittpunkten dieser Strahlenpaare bestehende 
„optische Modell‘ stellt sodann eine ähnliche Verkleinerung 
des aufgenommenen Geländes dar. 

Um ein Maß für die Abweichung zweier zugeordneter 
Zielstrahlen 8’, 8’’ von ihrer Schnittlage zu erhalten, hat 
man den Begriff der y-Parallaxe eingeführt. Zur näheren 
Erklärung sei ein rechtwinkliges Koordinatensystem an- 
genommen, und zwar — in Übereinstimmung mit dem 
Normblatt ‚Begriffe, Benennungen und Formelgrößen in 
der Photogrammetrie (Bildmessung)‘‘, DIN, 3035 Entwurf 
Juni 1951 — ein ‚‚Rechtssystem‘‘, dessen Ursprung im lin- 
ken Aufnahmezentrum 0’ liegt, dessen positive z-Achse das 
rechte Zentrum 0” enthält und dessen positive 2-Achse nach 
dem Zenith gerichtet ist (Bild1). Die Strahlen 8’, s’’ haben 
im allgemeinen eine einzige zur y-Achse parallelegemeinsame 
Bild Transversale t, die 8’, 8’ in zwei Punkten T', T’ schneidet. 

Die Länge der Strecke T’’'T’ heißt sodann y-Parallaze p, der 
beiden Strahlen 8’, 8’. Sie ist-positiv oder negativ, je nachdem die Richtung T’ — T’ mit der 
positiven y-Richtung übereinstimmt oder nicht übereinstimmt. Insbesondere kann p, = 0 
als Kennzeichen dafür gelten, daß s’ und s’’ einander im Raume schneiden; sie liegen dann 
gemeinsam in einer (durch 0’ und 0” gehenden) Kernebene. 

Wir denken uns vorerst den Orientierungsvorgang umgekehrt ausgeführt, nämlich beide 
Bündel, ausgehend von ihrer exakten gegenseitigen Orientierung, irgendwelchen Lagenände- 
rungen unterworfen, und fragen nach der dabei hervorgerufenen y-Parallaxe p, zweier zuge- 
ordneten Zielstrahlen 8’, 8”. Jede Verlagerung eines der Bündel läßt sich bekanntlich aus ein- 
fachen Grundbewegungen zusammensetzen, nämlich aus Verdrehungen um Achsen, die mit den 


*) Mitteilung aus dem Bundesamt für Eich- und Vermessungswesen in Wien. 

1) Gekürzte Wiedergabe von Vorträgen, die der Verfasser im Februar 1 
abgehalten hat. — Das gleiche Thema behandelte der Verfasser am 1. X 
für Photogrammetrie an der Technischen Hochschule München. 


953 im oben genannten Bundesamt 
II. 1953 in einem Vortrag im Institut 
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ET EEK Be 
Koordinatenachsen zusammenfallen oder dazu parallel sind, und aus darauffolgenden Ver- 
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schiebungen in den Richtungen dieser Achsen. Wir bezeichnen wie üblich die Drehungen um 


Achsen, die zur &-, y- oder z-Richtung parallel sind, in dieser Reihenfolge als Querneigungen 
(Kippungen), Längsneigungen (Verschwenkungen), bzw. als Verkantungen. Sind die zuge- 


hörigen Drehungswinkel w', &’ bzw. g', @’ bzw. x’, x'’ gegeben, ebenso die Verschiebungs- 
strecken in den Achsenrichtungen, so ist damit eine bestimmte gegenseitige Verlagerung der 


beiden Bündel festgelegt. DR : 
Zwei zugeordnete Zielstrahlen s’, 8’, die sich in der Ausgangslage in einem Punkt P des 


optischen Modells schneiden, erhalten nach der Bündelverlagerung eine y-Parallaxe 7p,, die 
sich unschwer als Funktion F der Koordinaten &, y, 2 dieses Punktes, der Modellbasis b und 
der gegebenen Größen w, - - - darstellen läßt. Für genügend kleine Werte dw, do’, dp’, dp”, 
dx’, dx’, dbz, dbz, db,, db,, db;, db; der 12 unabhängigen Veränderlichen dieser Funktion, wie 


solche meistens in Betracht zu ziehen sind, kann man F nach Taylor entwickeln und die Glieder 


zweiter und höherer Ordnung weglassen. Darnach ergibt sich als erzeugte y-Parallaxe: 
2 2 = 
dp, = ee (do — dw”’)-+ = dp! — si BY dp’ + dx’ | 
4 3 s ee, 
— (© 5) dx” + db, — db, — (db; — db; | 


Der Gültigkeitsbereich dieser Gleichung kann nach R. Finsterwalder?) für 
normal geöffnete Aufnahmen zwischen + 0,5® bis + 12, für Weitwinkelaufnahmen zwischen 
+t1® bis +28 der Bündeldrehungen angenommen werden, wobei die Grenzen der Bündelver- 
schiebungen auf entsprechende Drehungen zurückführbar sind. 
In der Praxis wird die gegenseitige Orientierung auf zwei verschiedene Arten herbeige- 
führt, nämlich entweder mit Festhalten der beiden Bündelscheitel im Raume oder so, daß eines 
der Bündel in Ruhe bleibt und nur das andere eine geeignete Verlagerung erfährt. Der erst- 
genannte Vorgang kommt zur Durchführung, wenn ein ‚„‚unabhängiges‘‘, d.h. ein einzelnes 
Bildpaar vorliegt, hingegen ist das zweite Verfahren, der sogenannte ‚‚Folgebildanschluß‘‘, 
wiederholt anzuwenden, wenn (wie bei der Aerotriangulation) die gegenseitige Orientierung 
aller aufeinander folgenden Bilder eines Bildstreifens herzustellen ist. 

Für diese Sonderfälle hat man in (1) anzunehmen: 
A) Bei „unabhängigen Bildpaaren‘“: 
do — do" —=do, dy,—=0, dd, =0, dd, =0, du, =0; 


“ B) beim ‚‚Folgebildanschluß‘‘, wenn das rechte Bündel festgehalten wird: 
ao = 0, net, te > 0, db, elchbzw. 


C) wenn das linke Bündel fest bleibt: 
de et dt det aueh 
In jedem dieser Fälle ist die Parallaxe dp, gemäß (1) als homogene lineare Funktion von fünf 
Größen dw, dp, : -- dargestellt, welche die Bündelverlagerung bestimmen. Die relative Quer- 


neigung dw’ — dw’’ der beiden Bündel wird der Einfachheit halber meist nur mit dw bezeichnet, 
wobei im Falle B) dw = dw’ bzw. für 0) do = — dw’’ zu setzen ist. 


Schreibt man in der allgemein gültigen Formel (1) zur Vereinfachung 


2.d#' — (2 — b) da’ + db, — db, = dy \ (2) 
— dp’ + (x — b) dp” + db, — db, = de Bl: 
so geht (1) über in 
2 2 
I ’—* dw Eu re 


Damit sind die von den einzelnen Grundbewegungen der Bündel herrührenden Parallaxen- 
anteile in einer Weise zusammengefaßt, daß hieraus umgekehrt die Gl. (1) leicht plausibel zu 
machen wäre. 
Um die letzte Gleichung geometrisch zu deuten, dividieren wir sie durch dw = dw — dw” 
(7 0) und setzen die dabei zu Tage tretenden Strecken (von nicht beschränkter Länge): 
dp, dz day 


us -n =: 5 TEE ne (AS 


®\R.Finsterwalder, Photogrammetrie, 2. Auflage, Berlin 1952, S. 156 


_ Davon haben Y, Z (ebenso wie dy, de in € 


betrachten nun Y, Zals y- 


a, \ Ä 

ausschlaggebender Bedeutung, aß 

schem Wege rasch gefunden werden 

Um dies darzulegen, geben wi 

(4) die Form: R - he 
| ET Ye 

Dies besagt offenbar, daß innerhalb jeder Ebene 


'orerst der Gl. 


mit Cden Punkt mit den Koordinaten Y, Z— 8, so steht die 
Gerade n—P-—( nach (5) immer normal auf dem Verbin- 
dungsstrahl vvon P mit dem Punkt B (z;, 0,0) auf der Basis. 
Dabei ist © der Gegenpunkt von B auf dem Kreis e. Da die 
Punkte 6(Y,Z— S)und 6(Y, Z)aufeiner Parallelen zurz-Achse 


ibn | liegen und den Abstand S haben, kann man sogleich feststellen : 


 Hilfssatz: Gegeben sei eine beliebige, jedoch in gewissen engen Grenzen gehaltene Bündelver- 


lagerung; wird durch sie bei einem Punkt P des optischen Modells, der die Koordinaten &,, Y, 4 


besitzt, die (genügend kleine) y- Parallaxe dp, hervorgerufen und bezeichnet dw die an dieser Ver- 
lagerung beteiligte Querneigung, so ergibt sich stets eine durch den Grundpunkt & der Ebene x = =, 
gehende Gerade q, wenn man im Punkt P die Normalen zur Kernebene y :z —= konst. = , : 2; (in 
der P liegt) errichtet und n in der z- Richtung durch die Strecke 8 — dp,/dw verschiebt (siehe Bild 2, 
wo Z und S negatives Vorzeichen haben). 

Da die Grundpunkte aller Normalebenen zur Basis wegen (2) und (4) auf einer Raumgeraden g liegen, er- 
füllen die in diesen Ebenen befindlichen Kreise € eine „‚orthogonale‘“ Regelfläche zweiten Grades. Auf dieser sind 
nämlich g und die Verbindungsgerade der beiden Zentren 0’, O’’ zwei zu reellenzyklischen Ebenen normale ‚‚Haupt- 
erzeugenden“. Die Fläche besteht aus allen Raumpunkten P, bei denen die betrachtete Bündelverlagerung dieselbe 
(genügend kleine) Parallaxe dp, hervorruft. Den verschiedenen Werten von dp, entsprechen ool gleichartige 
Flächen ‚konstanter Parallaxe‘“, die ein lineares Büsche] bilden. Damit läßt sich die räumliche Verteilung der 
Parallaxenwerte besonders gut veranschaulich>n. Die dp, = 0 entsprechende ‚„‚Grundfläche“ des Büschels umfa ßt 
insbesondere alle Raumpunkte, bei denen die Verlagerung keine Parallaxe erzeugt, so daß die gegenseitige Orien- 
tierung der beiden Aufnahmen einer solchen ‚gefährlichen‘ Fläche unbestimmt wird. Ähnliche Ergebnisse gelten 
auch für den allgemeinen Fall, daß die Drehwinkel und Schubstrecken der Bündel beliebige (unbeschränkte) 
Größen haben. Hierauf und auf eine Fülle weiterer Einzelheiten wird hier nicht eingegangen ®). Es sei nur er- 
wähnt, daß die Kenntnis dieser Zusammenhänge die Herleitung des in Rede stehenden graphischen Verfahrens 
wesentlich gefördert hat. s 

In der Praxis werden die beiden Zielstrahlbündel mit dem Einlegen der Bildplatten in das 
Auswertgerät wiederhergestellt und sogleich in eine Lage gebracht, die von der gesuchten 
gegenseitigen Orientierung nur mehr wenig abweicht. Hierauf mißt man an verschiedenen 
Stellen des optischen Modells die etwa noch vorhandenen y-Parallaxen 9; und hat sodann die 
beiden Zwelstrahlbündel (Projektoren) derart gegeneinander zu verlagern, daß alle beobachteten 
y-Parallaxen p; gleichzeitig zum Verschwinden kommen, oder anders ausgedrückt, daß bei allen 
betrachteten Modellpunkten jeweils die Parallaxe —p; erzeugt wird. 

Die diese Bündelverlagerung bestimmenden Orientierungsgrößen do’, do’, - - - genügen 
dem aus Gl. (1) folgendem System linearer Gleichungen, wenn dort nacheinander die Koordi- 
naten %, %, 2; der ausgewählten Modellpunkte samt den zugehörigen Parallaxen — p; einge- 
führt werden. 

Da die gegenseitige Orientierung sowohl beim ‚‚Folgebildanschluß‘‘ wie auch bei ‚‚unab- 
hängigen Bildpaaren‘“ durch fünf Grundbewegungen hergestellt wird (s. weiter oben) sind in 
beiden Fällen mindestens fünf Wertverbindungen %, Yi 2 Pi (£ 25) als bekannt vorauszu- 
setzen. Daraus ergeben sich im allgemeinen die fünf Orientierungsgrößen do, dp, --- durch 
Auflösen des entsprechenden Gleichungssystems. 

Anstatt dieses umständlichen numerischen Verfahrens kann die Hauptaufgabe auf Grund 
des erwähnten Hilfssatzes auch auf folgende Aufgabe der Liniengeometrie zurückgeführt wer- 
den *): „Gegeben sind fünf paarweise ähnliche Parallelstrahlbüschel, die in fünf untereinander 

®) Der allgemeinen Klarstellung dieses Fragenkreises hat J. Kra me seine Reihe von Veröffentlichungen 
gewidmet, auf die u. a. in den weiter unten zitierten Arbeiten hingewiesen wird. 


”) J. Krames, Gegenseitige Orientierung von Luftaufnahmen mittels liniengeometrischer Konstruk- 
tionen. Anzeiger d. Österr. Akad.d. Wiss., math.-nat. Kl. 86 (1949), S. 128—135. ‚ TE en 


bei allen Punkten P eines Kreises e_ dieselbe (kleine) Par Br “ 
dp, = 8 (do — do’) erzeugt wird (s. Bild 2). Denn bezeichnen 
wir mit 9, z die Koordinaten eines solchen Punktes P, ferner 


5. 
5 

eis Fr 
ud Ku 


ler: Sa en zur Er Die L 
praktische Anwendungen zu kompli 


vereinfachen, hat. man zoch die Ba pain Ei in mr % 


nzunehmen. 3 


- Wir zeigen vorerst, wie der Grundpunkt & jeder zur Basis normalen Ebene Se konst.=#;) 
auf einfachem graphischen Wege gefunden wird, sobald bei drei in dieser Ebene enthaltenen 


_ Modellpunkten P,, P,, P,die y-Parallaxen p,, P, ps gemessen sind. Zu diesem Zweck zeichnen 


wir die Verbindungsgeraden V,, V2, V, des Basispunktes B («,, 0,0) mit P,, P,, P, und errichten 


‚in diesen Punkten die zu V,, Vz bzw. v, normalen Geraden n,, n,, n, (Bild 3). Solange die zum 


Wegschaffen der Parallaxen p,, Ps, pz erforderliche Querneigung dw = dw’ — dw'’ nicht ge- 
funden ist, sind auch die (nach unserem Hilfssatz bestimmten) Der sppl DuIE EI cn en die 


jetzt zu setzen sind, zunächst unbekannt. Es gilt aber: 


dev 
SR 83 = P1: Pa: P3- 


Wir wählen nun einen positiven oder negativen Faktor F und verschieben die Normalen 
n,,N;, n, in der z-Richtung durch die Strecken K, = Fp,, K,= Fp,, K;, = Fp,, die den ge- 
gebenen Parallaxen, also auch den Strecken $; proportional sind. Die drei Normalen n,,n;,n, 
und ihre neuen Lagen g,, 9,, 4, nach den Verschiebungen bilden sodann zwev ähnliche und ähn- 
lich gelegene Dreiecke, und ihr Ähnlichkeitszentrum deckt sich bereits mit dem gesuchten Grund- 
punkt @. Denn nur an dieser Stelle ist der Hilfssatz für alle drei Punkte P; gleichzeitig erfüllt >). 


Z 


Bild 3 Bild 4 


Ist aber der Grundpunkt @ gefunden,.dann sind auch die in der 2-Richtung gemessenen 
Strecken S,; zwischen @ und den Normalen n, gegeben und die gesuchte Querneigung beträgt: 


! 77 Da ner Pass 
dev dv do st a 


Sorgt man dafür, daß in der Ausgangslage eine der drei gemessenen Parallaxen, etwa 
?,, = Oerscheint, dann liegt der Grundpunkt 6 auf der Normalen n, — q, (Bild 4). 


5) Aus diesen Überlegungen folgt zugleich, daß der Grundpunkt @ bereits eindeutig bestimmt ist, wenn bloß 
die Verhältnisse der drei Parallaxen p; bekannt sind. @ deckt sich nämlich mit jenem Punkt, dessen (mit den Fak- 
toren 1 :sin X zn; multiplizierten) „homogenen Koordinaten“ bezüglich des Grunddreiseits D,, D,, N, den ge- 
gebenen Parallaxen proportional sind. Vgl. hierzu J. Krames, Zur Geometrie der gegenseitigen Einpassung 
von Luftaufnahmen. Sitzgsber. Österr. Akad. d. Wiss., math.-nat., IIa, 160 (1951), S. 113—128. 
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Unter der Voraussetzung mathematisch exakter Strahlbündel ergibt sich aus jedem 
anderen Querschnitt & = konst., wenn darin drei Modellpunkte samt Parallaxen gegeben sind, 
wie oben gemäß (6) derselbe Wert von dw. Um daran anschließend auch die übrigen Orientie- 
rungsgrößen herauszufinden, sucht man insbesondere die Grundpunkte @’ und &” der Ebenen 
x —=0und = b. Aus den (2) und (4)entnommenen Koordinaten dieser Punkte: 

bdp’ — db; + db, ddr” + dby — dby, 


v zZ’ 


do ; dv - (7) 
yı _ bay — db + dbr gu _ bin’ + dby — dby | Ber. 
= dev : de 
folgen sofort die einfachen Formeln: 
A) für „unabhängige Bildpaare‘‘: 
[2 5 {£ Zn Be r. BL , F 
do = do, de = do, dx = do, dx -7dn. . (8); 
B) beim ‚‚Folgebildanschluß‘‘, wenn der rechte Projektor fest bleibt: 
dp - I a0, an - IE dr, dy- Lido, du=— Yda (9; 
C) beim ‚‚Folgebildanschluß‘“‘ mit festgehaltenem linken Projektor: 
dp" - IH anr, an -2—E dar, dd; = 2 de", diy=— Y" do" (10). 


Durch die beschriebenen graphischen Konstruktionen und die Gl. (6), (8), (9), (10) ist 
mithin die Hauptaufgabe für die erwähnte Voraussetzung gelöst ®). 

Hierzu seinoch bemerkt, daß die y-Parallaxen meistens in sechs charakteristischen Punkten 1, 2, 3, 4, 5, 6 
beobachtet werden, von denen 1, 3, 5 dem linken Querschnitt x = 0 und die übrigen 2, 4, 6 dem Querschnitt « =b 
angehören. Danach sind die Grundpunkte @’ und G’’ dieser Ebenen unmittelbar nach Bild 3 oder 4 zu finden. Zur 
weiteren Vereinfachung wählt man 1 und 2in den Schnittpunkten der durch die Zentren 0, 0°’ gelegten (z-par- 


Bild 5 


allelen) Lotstrahlen mit dem optischen Geländemodell (y,= 0, y; = 0), ferner die übrigen Punkte derart, daß die 
Verbindungsstrahlen v;, v; und V,, v, von 3, 5 mit 0° bzw. von 4, 6 mit 0” gegen die z-Richtung unter gleichen 
Winkeln + o geneigt sind (vgl. Bild 1). Zur genauen Festlegung solcher Punkte 3, 4, 5, 6 bedient man sich einer 
durchsichtigen Maske, auf der punktierte Kreisbogen gezeichnet sind (s. Bild 5). Diese wird abwechselnd auf den 
beiden (im Auswertgerät fixierten) Bildplatten in eine solche Lage gerückt, daß der gemeinsame Mittelpunkt M der 
Kreise jedesmal auf dem Lotstrahl liegt. Die Modellpunkte 8, 4 und 5, 6 sind sodann von den Betrachtungsoku- 
laren aus auf einem der Kreise einzustellen. Aus seinem Radius a und der Brennweite f folgt sogleich tgo =aj/f. 


' Wird überdies von einer Anfangslage ausgegangen, bei der die Parallaxen bei 1 und 2 be- 
reits eliminiert erscheinen, dann können die Grundpunkte @ und @” nach Bild 4 besonders. 


rasch konstruiert werden. Man legt dabei zweckmäßig die beiden Querschnitte &—=0 und 
& = b.derart übereinander, daß die Strahlen v,;, v4sowie Y,, vg zur Deckung kommen (s. Bild 6). 


6), J.Krames, Über ein graphisches Verfahren zum gegenseitigen Ein assen Luftaufnah 
Zeitschr. f. Verm. 37 (1949), 8. 1329. = : z nz 


se 13 _Parallaxen Dir nn si e 
In einem ‚solchen Falle en vor allem die beiden Werte 
ng überein, die sich — wie oben beschrieben — aus je drei Par- 
ınge in den Querschnitten & = 0 und 2—=bergeben. Dies trifft insbesondere 
el beiden Zielstrahlbündel exakt wiederhergestellt sind und alle Parallaxen allein von 
= enseitigen Lage der Bündel herrühren. 
Eine derartige Parallaxenverteilung liegt indessen nur sehr selten vor. Vielmehr sind die 
beiden Zielstrahlbündel vorwiegend mit solchen Beobachtungs- und Instrumentenfehlern be- 
haftet, daß dw, # dw,, ist. Danach liegt mit der Hauptaufgabe zugleich eine Aufgabe der 
Ausgleichsrechnung vor, d.h. die ursprüngliche gegenseitige Orientierung der Bündel kann 
- — überhaupt nicht exakt wiederhergestellt werden, sondern nur die „‚wahrscheinlichste‘‘i im Sinne 
5 des Prinzipes der kleinsten Quadrate von C.F. Gau. x 
Um diese Orientierung herauszufinden, können die beschriebenen graphischen Opera- Ur 
tionen unverändert beibehalten bleiben, womit der größte Teil der sonst erforderlichen Rechen- i 
u "arbeit umgangen wird. Hingegen dürfen die Gl. (6), (8), (9) und (10) jetzt nicht mehr ange- 
% wendet werden, sie sind vielmehr wie folgt umzuformen ®). 
a) Als Größe dw (bzw. dw’ oder dw’) wird in allen Gliedern von (8), (9), (10) durchgehend 
do, (oder auch dw,7) eingesetzt; überdies erhält J 


b) jede dieser Formeln ein additives Zusatzglied Be. 


5 (PER) (de De Re 


angehängt, Darin bedeuten J und X zwei Konstante, von denen die eine, etwa J, derart ange- 
nommen sei, daß aus den Gleichungen alle Produkte von dw, mit Y’” oder Z” (bzw. von don 
mit Y’ oder Z’) verschwinden. Wie sich alsbald zeigen wird, verfolgen diese Maßnahmen den : 
Zweck, jede der Orientierungsgrößen als lineare homogene Funktion der Parallaxen p, darzu- 
stellen. Wird beispielsweise in der nach a) und b) ergänzten Formel (9) für dx’ 


iz 


de’ = . dor +(J + K) (do; — don) der Faktor J = — - 4 

gewählt, so ergibt sich: & 

de’ = Z" Ber zZ’ = ER (dor don) 2 er (gi) . E 

R ce) Anstatt Gl. (6) setzen wir sinngemäß: 7: 
E ea ee IR OR k 


Alle Glieder (11) fallen weg, wenn Aw = dwı — do, = 0 ist, also wenn ‚‚geometrisch 
mögliche‘‘ Parallaxen vorliegen. 

d) Schließlich wird der in (11) enthaltene Faktor K für jede gesuchte Orientierungsgröße 
derart festgelegt, daß ihr mittlerer Fehler zu einem Minimum wird. Zu diesem Zweck formen 
wir (3) gemäß (4) um auf: 


yYınZtan ent We9,5 oder —2,4,6). zesart) 


. Darin gibt dw = dw’ — dw’ jetzt — im Gegensatz zu (1), (3) — jene relative Querneigung der 
beiden Bündel an, welche zum Wegschaffen der in den drei Punkten P; vorhandenen y-Par- 
allaxen p; erforderlich ist (vgl. auch Gl. 6). Setzt man in (12) der Reihe nach die Koordinaten 
%Y;, 2; der in x = 0 gelegenen Modellpunkte 1, 3, 5 und die zugehörigen Parallaxenwerte 2, ein, 
so erhält man drei lineare Gleichungen für den reziproken Wert der Querneigung dw, und die 
Koordinaten Y’, Z’ des Grundpunktes @'. 


Mit der vereinfachten Bezeichnung der Determinanten: 
yi+aia ma =UÜ, vv Yti »nal=TV 
| yı  pal=M, M% y+al=N, 
folgt daher aus (12): 


vr FrrUS 7 =E dt do, = =, Pe EI ee (13). 


?) Man erkennt diesauch daran, daßim System der sechs entsprechenden linearen Gleichungen (s. oben) nur 
fünf linear unabhängige vorhanden sind. 

8) J. Krames, Erweiterung des graphischen Einpassens von Luftaufnahmen auf den Fall eines vorliegen- 
den Widerspruches. Schweiz. Zeitschr..f. Verm. u. K. 49 (1951), S. 293—299, 307 —317. 
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‚ergibt sich z. B. aus (9*) oki 3 Er ws. 

# sell N ee ie SER “ s Alt x 
2.0, We N 
Mithin sind tatsächlich alle Orientierungsgrößen gemäß a), b) und e) als lineare h mogene 
Funktionen der sechs y-Parallaxen p; gegeben. er a Or 

Gleiches trifft bekanntlich auch zu, wenn diese Größen aus den ‚‚vermittelnden Beob- FE 
achtungen“ 9; mit Hilfe von ‚‚Normalgleichungen‘‘ berechnet werden. Um beide Darstellungen 
zur Übereinstimmung zu bringen, ist nur noch der Faktor K in geeigneter Weise festzulegen. 
‘ Denn für Aw = do, — dor = 0 soll ja jede Gleichung von der Form (15) in die entsprechende 

in (8), (9) oder (10) übergehen, deren rechte Seite wegen (13) und (14) ebenfalls eine lineare 
homogene Funktion der Parallaxen », darstellt. Diese Forderung wäre mit keinem anderen » 
Be: Ansatz als (9*) erfüllbar. i 
er: Zur Berechnung der Faktoren K schreiben wir die gemäß a), b), e) umgeformten Gin. (8) 
(9), (10) und (6*) einheitlich in der Form : “ 


el + Kup. Pik- LOS 


Hierin ist da die betreffende Orientierungsgröße, während s;, tx Konstante bedeuten. Aus (16) 
folgt als ‚, Gewichtskehrwert‘‘ von da: u 


oda \? ’ : 
On = 2 (on) = I (+ +Ku? = (x +2K st, + K?ti). 


Um sein Minimum zu erhalten, setzen wir 


d ' 
Des 2 2 gl + 2K 2-0 
und beachten, daß (für reelle Modellpunkte) stets gilt: 
A) Sa =>0 
dK? SE 
Also ist 
' Kz— N st: Dt ee a en a 


jener eindeutig bestimmte Wert des Faktors K, für den der mittlere Fehler m, = u aa von 
da am kleinsten wird (dabei gibt u den mittleren Fehler der Parallaxenmessung an). Mit K ist 
daher der wahrscheinlichste Wert jeder Orientierungsgröße eindeutig bestimmt. Auf ähnliche 
Weise — ohne Auflösen von Normalgleichungen — wurden diese Werte anscheinend noch 
nirgends berechnet. Wie aus (12) bis (16) leicht ersichtlich ist, sind alle Faktoren X allein durch 
die Raumlage der Punkte 1,2, 3, 4, 5, 6, also durch die Form des aufgenommenen Geländes ein- 
deutig bestimmt. 

Wir erteilen noch den sechs Modellpunkten 1, 2, 3, 4, 5, 6 die mit ER LER 
Yy:3a—- ya — —Y:2—= — Y:2,— tg o gekennzeichneten besonderen Lagen (vgl. wie- 
der Bild 1) und berechnen hierfür die oben verwendeten Determinanten. Nach dem Ein- 
führen der dabei erhaltenen Werte sowie der gemäß (17) berechneten Faktoren in die nach (15) 
erweiterten Gl. (6*), (8), (9), (10) ergeben sich sodann die gesuchten Endformeln für die Orien- 
tierungsgrößen. Dabei treten mehrere neue Größen hervor, die jedoch aus der Konstruktions- 
zeichnung (Bild 6) ebenfalls entnommen oder rasch abgeleitet werden können. Dies gilt vor 
allem für: 

tea +1=1/cos?o =e, 


Mr eu ren =l,, A . 


kg ey =U rg TER —=Us u+u=u, 


EN Dr le Li, a 


N EEE I A EN u ET EI EEE > 0) TE nn 
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desgleichen für die y-Koordinaten n’ und n’’ der Schnittpunkte N’ 
und N’ der Normalen n,, n, bzw. n,,n,; (s. Bild 6). Wir ermitteln 
ferner mit Hilfe eines Kreises vom Durchmesser 1 (Bild 7) die 
Faktoren : 

w2 w2 


> w?+ w'2 ’ L = Ww2+ ww’? 


Re) 


und führen folgende Abkürzungen ein: 


LE ne u a lea 
Das Gesamtergebnis lautet danach: 
A) Unabhängige Bildpaare: 
dv —= Rdo, + T do; , 
‚ _ Y" don + Rn’ Aw 


„ _.Y do, — Ty Aw 
dp 7 dp ; 
ne Z’ day + RE’ Aw „ a 2ZA0d, LEAD | ee 
Me ’ dx = 5 
b b J 
B) Folgebildanschluß, rechter Projektor fest: 
do’ = Rdw, + Tdwr, 
en a a A 
due a REO, = Z di TE Ai, 
. (20). 


‚_ 2’ don — 2’ dor +L Aw 

— z i 
C) Folgebildanschluß, linker Projektor fest: 
dw’ = Rdo, + T don: 


dx 


db, = — Y’ do, + Tn' Aw 


„ 0 d 1 2% d x A „ 7 
do I O1 z & +n 0} db, En di, AL RER Aw, ) 
Zr do“ rt do” +5 40 = £ ; . (21). 
dx — E 5 . i db, = — Y’ dv — Rn’ Aw | 
Für die Querneigungen gilt darin allgemein: 
a a Rn. 2. 


Obige Formeln bleiben im wesentlichen ungeändert, wenn in der Ausgangslage bei den 
Modellpunkten 1 und 2 keine Parallaxen mehr zu beobachten sind (vgl. oben S. 257, 258). 
Es fallen dann bloß die ersten Glieder in den Zählern von (22) weg und es gilt (vgl. Bild 4 
und 6): Z =z,Z”’=z,. Damit ist aber das auf allgemeine Parallaxenangaben erweiterte 
graphische Verfahren auf seine einfachste Form gebracht. 


Mit der Ableitung dieser Ergebnisse war die Entwicklung desin Rede stehenden Verfahrens keineswegs ab- 
geschlossen. In theoretischer wie auch in praktischer Hinsicht waren noch mannigfache Ergänzungen auszu- 
arbeiten. U. a. wurde klargestellt, inwieweit ein allfällig auftretender kleiner Winkel zwischen der Aufnahmebasis, 
d.i. die Verbindungsgerade der beiden Projektionszentren, und der &- oder z-Richtung des Maschinensystems zu 
berücksichtigen ist. Eine eingehende Untersuchung galt ferner jenen Fällen, bei denen die sechs Orientierungs- 
punkte nicht innerhalb von x = 0 und & = b gewählt werden können, sondern nur (zu je dreien) in zwei anderen 
Normalebenen zur Basis ®). Dies kommt mitunter vor, wenn das aufgenommene Gelände Wasserflächen oder 
Schattenpartien enthält oder wenn beide Aufnahmen sich zu wenig übergreifen. 

Um das Verfahren für den ständigen Gebrauch in der Praxis geeignet zu machen, mußte es schließlich den 
verwendeten Geräten tunlichst angepaßt und weitgehend schematisiert werden. Über seine Vorzüge und die bei 
seiner praktischen Handhabung gesammelten Erfahrungen wird bei anderer Gelegenheit berichtet. 


®») J. Krames, Graphisches Einpassen von. Luftaufnahmen bei beschränktem Gesichtsfeld samt An- 
wendung auf ein praktisches Beispiel. Österr. Zeitschr. f. Verm. 41 (1953), 2. Heft, S.1—9. 


Eingegangen am 13. Juni 1953. 
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U us Programme für Rechenautomaten durch Verm 
und zu vereinfachen, wird deren Struktur untersucht 


7. 


The general structure of programs for aulomalic compuling machinery is ö vestigat 
abbreviating and simplifying programs by eliminating redundancy of information. It 


can be rendered very versalile by decomposing it into elementary orders 
An example of an order code generated by such elementary orders is given. ion 
of methods for automatic control of changing processes is discussed in detail. en 
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$ 1. Einleitung 


Bei der Entwicklung programmgesteuerter Rechenmaschinen geht selbstverständlich 

das Bestreben dahin, die Programmierung im Regelfall kurz und übersichtlich ausführen zu. 
"können, ohne daß dabei die Anpassungsfähigkeit an verschiedenste Aufgaben im Rahmen eines 

abgesteckten Aufgabenbereichs beeinträchtigt wird. Übersichtlichkeit und Kürze vermindern 

die Gefahr von Fehlprogrammierungen und erleichtern dem Mathematiker die Bewältigung 

umfangreicher Aufgaben. Die konkrete Kürze des fertigen Programms wirkt sich überdies 

günstig auf Eingabezeit, Rechendauer und benötigte Speicherkapazität aus. , 

In der vorliegenden Arbeit wird nach einer Voruntersuchung über die allgemeine Struk- 
tur von Rechenprogrammen ($ 2) versucht, diesem Ziel der Verkürzung und übersichtlicheren 
Gestaltung der Programme durch geeignete Maßnahmen näher zu kommen, nämlich: 

a) durch Automatisierung gewisser strukturell bedingter Änderungsprozesse, die beim 
Zusammenbau (Änderungen erster und zweiter Art)und beim wiederholten Ablauf (Änderungen 
dritter Art) von Teilprogrammen auftreten ($ 3), 

b) durch zweckmäßige organisatorische Aufgliederung des Befehlscode in Befehlsele- 
mente ($ 4). 

Die Maßnahme a) befreit den Zusammenbau von Teilprogrammen vom Ballast umständ- 
licher Vorbereitungsprogramme. Die Maßnahme b) entlastet die Teilprogramme selbst, indem 
sie gewisse Mitteilungen an die Maschine überflüssig macht. 

Die Durchführung der vorgeschlagenen Maßnahmen wird in ihren Grundzügen an Bei- 
spielen erläutert ($5). In einem Anhang werden im einzelnen Möglichkeiten für die Befehls- 
steuerung zur Automatisierung der Änderungsprozesse besprochen. 


$ 2. Allgemeine Struktur der Rechenprogramme 


2.1. Teilprogramme undahre Gliederung 


Numerische Verfahren zerfallen im allgemeinen in eine Reihe einzelner Algorithmen. . 
Diese Eigenschaft überträgt sich naturgemäß auf die zugehörigen Rechenprogramme. Dem- 
entsprechend setzen wir fest: ‚‚Teilprogramm‘“‘ im Rahmen eines Rechenprogramms ist eine 
Befehlsfolge, die die Übersetzung eines bestimmten abgeschlossenen Algorithmus in die 
Sprache der Maschine darstellt‘ Damit ist das Teilprogramm eine neue, vorzugsweise rekursive 


A Arbeitsgruppe für elektronische Rechenanlagen der Technischen Hochschule München. 
) Die Arbeit wurde unterstützt durch ein Forschungsstipendium der Deutschen Forschungsgemeinschaft, 


u an nad 21 200; 


löse Den mathematischer en in einer ge die der üblichen völlig 


gleichberechtigt ist. 


. Das Teilprogramm weist eine topologische Struktur auf, die sich aus einer kleinen A 
von Elementen (einfache Strecke, Verzweigung von echt alternativem Charakter, Iterations- 
schleife und Induktionsschleife) ®) aufbaut und durch die bekannten Flußdiagramme ver- 


 anschaulicht wird. Diese topologische Struktur liefert die Verknüpfung der einzelnen arith- 


metischen Bestandteile des Teilprogramms und ergibt mit ihnen zusammen die Operatordar- 
stellung des Algorithmus. Sie findet ihren Ausdruck in Strukturbefehlen *)» die im 
allgemeinen im Teilprogramm selbst operieren und deshalb in ihm festliegen. Die einzige Aus- 
nahme stellen die Ausgangsbefehle dar, über deren Ziele man selbstverständlich willkürlich 
verfügen können muß, da sie die Verknüpfung mit anderen Teilprogrammen herstellen. 

Als Bestandteile des Algorithmus und damit des Teilprogramms selbst (nicht us Ele- 
‚mente seines Operationsbereichs) kommen zwei Arten von Zahlen vor: 

a) Rechengrößen, die den arithmetischen Bestandteilen des Algorithmus Ereelaheh 
sind, zum Beispiel feste Gewichte bei gewogenen Mitteln, wie etwa bei der Si impsonschen 
Integralformel. 

b) Parameter, die topologischen Strukturen zugeordnet sind, vornehmlich Durchgangs- 
zahlen für Schleifen. Für Induktionsschleifen sind sie gesondert vorzugeben und im allge- 
meinen beliebig variierbar. Bei Iterationsschleifen dagegen werden sie durch eine Rechengröße, 
die ebenfalls gesondert vorzugebende Genauigkeitsgrenze, ersetzt. 

Die Rolle der Variablen der mathematischen Formel wird von den Adressenteilen 
Datumadressen) der Befehle arithmetischer Programmbestandteile (meist arithmetische und 
Transfer-Befehle) übernommen, denn die zugehörigen Adresseninhalte, die Rechengrößen, sind 
der Operationsbereich des Algorithmus. 

Aus praktischen Erwägungen heraus wollen wir fordern, daß diese Rechensrößen nor- 
malerweise völlig unabhängig vom Teilprogramm gespeichert werden sollen und nicht etwa in 
eigens dafür vorgesehenen Leerstellen des Teilprogramms, das demnach ausschließlich die Be- 
fehlsfolge des abstrakten Algorithmus und gegebenenfalls Zahlen der oben beschriebenen beiden 
Arten enthält. Diese Regelung bringt gewisse Vorteile, zum Beispiel, wenn ein umfangreiches 
Gesamtprogramm für eine größere Anzahl von Rechengrößen-Sätzen durchgerechnet werden 
soll, für deren Auswechselung. Sie ist auch praktisch nicht zu umgehen, wenn Induktions- 
schleifen auftreten, die mehrere Sätze von Rechengrößen verarbeiten, wie etwa beim Skalar- 
produkt zweier Vektoren. 


22.Änderungen erster und zweiter Art 


Für die nachfolgende Untersuchung wollen wir uns ausdrücklich auf solche Änderungen 
eines Teilprogramms beschränken, die den Algorithmus selbst nicht stören. Eingriffe in die 
Befehlsfolge oder die Operationsteile einzelner Befehle, die im übrigen nur selten von Nutzen 
sein dürften, sollen also von der Betrachtung ausgeschlossen werden. Damit kommen über- 
haupt nur noch Änderungen von Adressen und von Parametern vom Charakter der Durch- 
gangszahlen in Frage, die sich formal ebenfalls als Adressen darstellen lassen. 

Über diese Parameter, über die Adressen der Rechengrößen sowie über die Datumadressen 
der Befehle, die aus dem Teilprogramm herausführen, kann frei verfügt werden. Sie sollen 
daher hier in Anlehnung an Goldstine und v. Neumann’), jedoch in etwas abge- 
ändertem Sinn, als „freie Variable‘ bezeichnet werden. Eine Sonderstellung nimmt 
unter ihnen der Eingang, das heißt der erste Befehl des Teilprogramms, ein, dessen Lage bei 
der Eingabe willkürlich gewählt werden kann, aber im Teilprogramm explizit überhaupt nicht 
auftritt. 

Die Strukturbefehle dagegen, die im Teilprogramm selbst operieren und damit im Teil- 
programm liegende Adressen aufweisen, können nicht mehr willkürlich verändert werden, da 
ihre Adressen relativ zum ersten Befehl des Teilprogramms festliegen. Sie sollen daher hier als 
„gebundene Variable‘ bezeichnet werden. 

An den gebundenen Variablen ist augenscheinlich nur eine einmalige Translation um die 
Differenz zwischen tatsächlicher und programmierter Adresse des Anfangsbefehls des Teil- 
programms vorzunehmen. Es handelt sich also um die Goldstine-v.Neumannschen 
„Änderungen erster Art“ (GN 11.111.12.4—6). 


3) Die beiden letzteren unterscheiden sich im wesentlichen dadurch, daß bei der Iterationsschleife die Zahl 
der Durchgänge von: vorausgehenden Ergebnissen bestimmt wird, während sie bei der Induktionsschleife als 
Durchgangszahl explizit angegeben werden muß. 

4) Dies sind vornehmlich nichtaritbmetische Befehle wie Sprung und Adressensubstituticn. 

5) H. H. Goldstine — J. von Neumann, Planning and Coding of Problems for an Electronic 
Computing Instrument, Princeton 1948, im folgenden zitiert als GN unter Angabe von Teil und Band. 


feste Zahl, die dem Abstand der Komponenten im Speicher entspr 


(GN loc. eit.) unterliegen. 
Bez er 


Die vorangegangene Diskussion hat die Maı 
programmen noch nicht ausgeschöpft. Sie betrifft 


nacheinander auf. 


gebenen Endwert heraufgezählt werden, oder andere Maßnahmen, die das Abbrechen der 


Induktion veranlassen. 

Änderungen an Indexvariablen, die also vom Teilprogramm selbst nach jedem einzelnen 
Durchgang durch eine Induktionsschleife ausgeführt werden müssen, wollen wir „Ande- 
rungen dritter Art‘ nennen. TR 


24. Zusammenfassung 
Damit sind.nun drei durchaus verschiedenartige Gruppen von systematischen Ände- 
rungen festgestellt: i 
Die Änderungen erster Art sind einmalige Translationen, die ausgeführt werden 
durch Addition einer festen Adresse (zweckmäßig derjenigen des ersten Befehls des Teilpro- 
gramms im Speicher) auf die ursprünglich programmierten Adressen, die gebundene Variable 
darstellen. 


Die Änderungen zweiter Art sind Substitutionen frei wählbarer Adressen in die irrele- 


vant besetzten Adressenstellen, die freien Variablen entsprechen. Sie sind normalerweise bei 
jedem Aufruf des Teilprogramms vorzunehmen. ; 

Die Änderungen dritter Art zählen die Komponenten von Indexvariablen ab. Sie be- 
wirken dies im allgemeinen durch wiederholte Addition fester Zahlen auf bestimmte, meist 
wohl durch Änderungen zweiter Art fixierte Adressen. 


$ 3. Prinzipien der Automatisierung der Änderungsprozesse 


| 31. Vorbemerkung 
Wegen der Verschiedenheit der aufgeführten Änderungsprozesse ist eine automatische 
Durchführung nach einem einheitlichen Verfahren nicht zweckmäßig. Überdies gibt es im 
Ablauf eines Programms gerade dreierlei verschiedene Zeitpunkte, in denen die Maschine 
automatisch und ohne wesentliche Störung des Programmablaufs die Notwendigkeit der 
entsprechenden Änderungen feststellen und ihre sofortige Ausführung einleiten kann. 


32. Änderungen erster Art 


Für Änderungen erster Art ist dieser günstigste Zeitpunkt gekommen, wenn der erste 
Befehl des Teilprogramms vom Eingabemedium abgelesen und in den Speicher transferiert 
wird. Der Befehlszähler, der die sukzessive Einspeisung der Befehle in den Speicher regelt, 
enthält dann gerade die Translationsadresse, die zu sämtlichen gebundenen Variablen des sich 
anschließenden Teilprogramms zu addieren ist. Auf ein (nicht zu speicherndes) Anfangszeichen 
im Eingabemedium wird die Translationsadresse festgehalten. Sie wird auf alle nachfolgenden 
Adressen addiert, die durch ein zweites (ebenfalls nicht zu speicherndes) Zeichen als gebundene 
Variable markiert sind. Dies ist im wesentlichen das bei der EDSAG verwandte Verfahren. 
Unzweckmäßig dürfte es sein, die Änderungen erster Art nach erfolgter Speicherung der Teil- 
programme auszuführen. 


33. Änderungen zweiter Art 


as Für Änderungen zweiter Art kommt der günstigste Zeitpunkt der Ausführung, wenn der- 
jenige Sprungbefehl des Gesamtprogramms im Befehlssteuerwerk die Kontrolle übernimmt, 
der das fragliche Teilprogramm aufruft. Er kann und muß die neu einzuführenden Adressen 


Parallel laufen Änderungen der Induktionsindizes, die schrittweise bis zu dem vorge- 2% 
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sozusagen mitbringen. Der Gedanke liegt nahe, auch hier diejenigen Befehle in Bernie 


Weise zu kennzeichnen, die freie Variable enthalten. Sämtliche Befehle des Teilprogramms 
sind der Reihe nach auf dieses Kennzeichen hin zu prüfen und ihre Datumadressen gegebenen- 
falls sofort zu ersetzen. Die Prüfung betrifft jetzt Befehlsworte, die sich im Speicher befinden. 
Es ist daher zweckmäßig, sie im Rechenwerk vorzunehmen, das auch die Änderungen, nämlich 
Adressensubstitutionen, auszuführen hat. 

‚Die zu substituierenden Adressen müssen in irgendwelchen Speicherzellen greifbar sein. 
Die Adressen dieser Zellen wiederum müssen dem Leitwerk in dem Augenblick, in dem die Not- 
wendigkeit der einzelnen Anderung festgestellt wird, bekannt sein. Tatsächlich steht nun dem 
Leitwerk in dem Augenblick, in dem das Teilprogramm aufgerufen wird, nicht nur die Adresse 
des Anfangs des Teilprogramms (als Datumadresse im aufrufenden Sprungbefehl) zur Ver- 
fügung, sondern auch noch die Adresse der Speicherzelle, aus der der Sprungbefehl abgelesen 
wurde (im Befehlszähler). Ordnet man nun die der Reihe nach zu substituierenden Adressen 
in den unmittelbar an die zuletzt genannte anschließenden Speicherzellen an, so hat das Be- 
fehlssteuerwerk nur noch dafür zu sorgen, daß diese Sequenz an die gekennzeichneten Stellen 
des Teilprogramms einsortiert wird, während dieses Befehl für Befehl das Rechenwerk durch- 
läuft. Durch ein Schlußzeichen im letzten Befehl des Teilprogramms wird der Prüfvorgang 
abgebrochen und sofort die Ausführung des eben richtiggestellten Teilprogramms eingeleitet. 


3.4.Anderungen dritter Art 

Änderungen dritter Art schließlich sollen zu dem Zeitpunkt ausgeführt werden, in dem 
der zu ändernde Befehl im Befehlssteuerwerk eintrifft. In diesem muß dann auch die auto- 
matische Korrektur vorgenommen werden, zumal da das Rechenwerk jetzt fortwährend für die 
eigentliche Zahlenrechnung freigehalten werden muß. Die Korrektur kann bewerkstelligt wer- 
den durch Angliederung eines kleinen Addierwerks an das Befehlssteuerwerk, dessen Akkumu- 
lator stets die aus dem Speicher abgelesene Adresse aufnimmt, während das Addendenregister 
die Abzählung der Schleifendurchgänge vornimmt und damit automatisch auf die zu Index- 
variablen zu addierende Zahl eingestellt ist. Die Indexvariablen selbst sind durch ein be- 
sonderes Zeichen im Operationsteil des Befehls gekennzeichnet, durch das die Addition aus- 
gelöst wird. Unmittelbar anschließend kann die nunmehr richtig bezeichnete Speicherzelle 
aufgesucht werden. 

Soweit uns bekannt ist, dient das ‚‚Schrittregister‘‘ des Göttinger Rechenmaschinenpro- 
jekts G2 ähnlichen Zwecken. Die konstruktiv sehr ähnlichen sog. Indexregister, wie sie zum 
Beispiel in dem Bericht von H. Rutishauser über automatische Rechenplanfertigung 
(Nr. 3 der Mitteilungen aus dem Institut für angewandte Mathematik an der ETH Zürich, 
Basel 1952) beschrieben sind, stellen dagegen einen durch die Trennung von Befehls- und 
Zahlenspeicher bedingten Ersatz für die an sich erheblich einfachere Adressensubstitution dar. 
Sie sind nicht für die automatische Abzählung der Indizes eingerichtet. 

Es muß durchaus als möglich angesehen werden, daß Programme mit komplizierter topo- 
logischer Struktur mit den angegebenen Maßnahmen nicht mehr automatisch behandelt wer- 
den können. In solchen Fällen wird man auf explizite Programmierung der Adressensubstitu- 
tionen zurückgreifen. In der Regel jedoch wird die angestrebte Verkürzung der Programme 
und Vereinfachung der Programmierungsarbeit sicher zu erreichen sein. 


$4. Aufgliederung eines Befehlssystems in Elementarbefehle 
4, Elementarbefehle und VereineVon Elementarbefehlen 

Die Forderung der Übersichtlichkeit und Kürze erstreckt sich auf das ganze Programm, 
also über die Durchführung der bedeutsamen, aber speziellen organisatorischen Maßnahmen 
hinaus, die beim Zusammenbau von Teilprogrammen auftreten und von der topologischen 
Struktur herrühren, auf das Teilprogramm selbst. Sie betrifft dabei die Verkettung arithmeti- 
scher und organisatorischer Operationen, die die lokale Struktur der Befehlsfolge im Teilpro- 
gramm ausmacht, und läuft darauf hinaus, im Rahmen vernünftiger technischer Möglichkeiten 
das Teilprogramm mit möglichst wenigen Digits zu codieren ®). 

Da die Adressenstellen einen merklichen Anteil der insgesamt benötigten Digits aus- 
machen, hat man anzustreben, sämtliche zwischen zwei Adressen auszuführenden Operationen 
im eigentlichen Befehlsteil eines Befehlswortes auszudrücken ?). 

Das bedeutet, daß im eigentlichen Befehlsteil eines (Ein-Adreß)-Befehls Kombinationen 
gewisser Einzelbefehle vorkommen. Als „Elementar befehl‘ einer solchen Aufgliede- 


6) In der Sprache der Informationstheorie: von „redundancy‘‘ möglichst zu befreien. 3 
?) Wir denken dabei an ein Ein-Adreß-System. Bekanntlich zeigen Mehr-Adreß-Systeme diesem gegenüber 
eine Verschwendung von Digits. 


; A AC oder een ZEN 
b) bilde den Betrag der Zahl in Zeit n Et 1 
c) addiere oder subtrahiere die so bel 


2... Um zu einer gewissen , Übersicht zu kommen, 
in ‚Vereinen‘ zusammenzufassen, so daß nur Mitgli 
wenigstens teilweise der direkten Kombination fähig si 
: glieder enthalten, die sinnvollerweise innerhalb eines e 
- müssen, praktisch vielleicht von seltenen Ausnahmen abgesehe e 

glied „keine Wirkung‘ als neutrales Element („—“) angehören. F & 

einen Verein von nur zwei Mitgliedern: „„Affirmative Aussage‘‘ und an Wirkung en 


42. Beispiele mei N 
Um ein Beispiel zu geben, seien einmal nur die Speziesoperationen aa Eine f 
günstige Aufgliederung ist die folgende °) un 1-7 
E ‘1. Verein > mn \ 2 j Be “ 


2. Verein {mp > MD 


Ka [MD] nach AC 
2 3. Verein $ lösche AC, add [MD] nach a0 
m. multipl [MD] mit [MR] ar 
5 divid [AC] durch [MD]. 


Hier sind sämtliche zwanzig Kombinationen sinnvoll. Da Sprungbefehle mit Spezies nicht 

r : kombinieren, könnte man sie in den dritten Verein mit aufnehmen. Die Kombinationsmög- 

Be =. lichkeiten mit dem ersten und zweiten Verein werden dann allerdings nicht ausgeschöpft. Dies 

. ist ein ganz allgemeiner Zug jeder derartigen Aufgliederung üblicher Befehlssysteme. Auch ist 
die Aufgliederung und Zusammenfassung zu Vereinen von Elementarbefehlen in gewissen 
Grenzen willkürlich und soll hier zunächst mehr als methodische Richtlinie denn als starre 
Vorschrift angesehen werden. Sie kann dazu dienen, die Menge und Anzahl aller im eigent- 
lichen Befehlsteil eines Befehlswortes vorkommenden kombinierten Befehle zu umreißen. Deren 
Anzahl überschreitet selbst den Umfang ausführlicher Befehlslisten, wie die des Raytheon 
Computer, beträchtlich. 


4.3. Organisation der Aufgliederung in Vereine von Riemen 2225 
befehlen 


Es ist:nun keineswegs erforderlich, jedem kombinierten Befehl wie üblich ein unab- 
hängiges Codezeichen zuzuweisen. Für eine Codierung des Befehls stehen, wenn man sich etwa 
an die tetradenweise Eingabe von (dezimal verschlüsselten) Ziffern anlehnt, ebenfalls Tetraden 
von Digits zur Verfügung. Die sechzehn Digitkombinationen mögen mit irgendwelchen passend 
gewählten Zeichen bezeichnet sein. Rein äußerlich treten also ebenfalls direkte Kombinationen 
von Code-Vereinen von je sechzehn Elementen auf und es liegt nahe, das Befehlssystem ent- 
sprechend in mehreren Code-Vereinen von je höchstens sechzehn Mitgliedern einschließlich des 
neutralen Elementes darzustellen (eventuell durch Zusammenfassung von mehreren Grund- 
vereinen). 

Über die technische Ausnützung einer noch vorhandenen weiteren Untergliederung wird 
noch zu sprechen sein (vgl. Ziff. 5.2). Für die Programmierung ist sie unerheblich. 

Die Zusammenfassung geschieht vorteilhaft nach allgemeinen, einprägsamen Gesichts- 
punkten. Man wird etwa das Gerüst der Spezies- und Sprungbefehle in einem Code-Verein 
zusammenfassen und andere Code-Vereine für die Bestimmung diskriminierender Merkmale 


®) Die Auflösung solcher Elementarbefehle in logistische Einzeloperationen der Schaltkreistechnik ist 
Sache des technischen Aufbaues des Befehlssteuerwerkes. 

®) Hier wie im folgenden bedeute: Eine eckige Klammer [.. .] ‚Inhalt von... “, m die Speicherzelle, die 
der Adressenteil des Befehlswortes angibt. 


$ 5. Beispiele und Erläuterungen 


5.1.Beispielfür einen Drei-Zeichen-Code Ei; 
gr Für eine verhältnismäßig einfache Maschine mit festem Komma möge etwa die innere Bu: 
- Gliederung des Befehlssystems in drei Code-Vereinen mit je höchstens sechzehn Mitgliedern Bl 5: 
(also durch Tetraden von Digits Jarstelibar) Ausdruck finden. 
Erstes Code-Zeichen _ Erster Verein 
7 wirkungslos (,, [MD] — MD“) 
— [MD] — MD 
[MD]| > MD 
.— |I[MD]| > MD 
[m] — MD hi 
— [m] —MD ; 
E g — [ml > MD | Bi; 
Zweites Code-Zeichen Zweiter Verein ; 
be > Sr) wirkungslos 
' -F add MD nach AC- 
CE lösche AO, add MD nach AC 
BE .. BAT . mult MD mit MR 
E divid A0durch MD _ | E 
E Sprung, wenn AC 0 : 
Sprung, wenn AÜ <O 
unbedingter Sprung 
q-malige Rechtsverschiebung in AC!P) 
q-malige Linksverschiebung in AC 
: Stop 


Drittes Code-Zeichen Dritter Verein 
nach MR nach AC nach m Zum Abdruck 
[MR] [4C] — — 
[AC] [MR] er Bi 
[MR] [A] e [40] 
[40] [MR] —_ IMR] - 
[MR] [Ac] [40] — 
[40] [MR] [MR] 
[MR] [40] [40] [4C] 
k [AC] [MR] [MR] [MR] 
Im zweiten Verein ist noch für bis zu fünf weiteren Elementarbefehlen Platz. Hier 
mögen Sonderbefehle eingefügt werden, wie Adressensubstitution und spezielle Befehle zur 


Konvertierung, die nicht in Kombinationen vorkommen und für die Diskussion hier entbehr- 
lich sind. j 


SU om. 


} nach Zelle m je 


N, 


a 1 ha nn a oe 
NO om om 


BETEN 


52.Vorläufiges über technische Durchführung 


Die Aufgliederung ist so vorgenommen, daß der erste und der dritte Verein aus direkten 
Kombinationen von je drei echten Alternativen bestehen. Hier liegt eine Untergliederung der 


10) Angezeigt etwa durch 2-1 im Adressenteil des Befehlswortes. 


2 Die Aufgliederung 
zelnen technischen E} 


00.) | +] Transferiere m na 
(2) | = RN 
(9) | «[Kompiementiere ) MD 


ee Spezies, Sprungbefehle usw. 


4 R (9) | Vertausche AC mi MR A% En N Ban 
E> (6) | Transferiere *) AC nach Zellem | 3. Verein. BE : = E:- 
| (7) | + Bringe AC zum Abdruck 5 | | ep EN RE 

{ Br 


«| Komplementiere *) MD, 


*) nunmehrigen. 


53. Diskussion spezieller wichtiger Kombinationen 
I Falls die technischen Elementarbefehle (1) und (6) in ein- und demselben Befehlswort 
Be auftreten, beziehen sie sich selbstverständlich auf die gleiche Adresse. In Verbindung mit 


5 Spezies-Befehlen ergibt sich damit die praktisch wichtige Möglichkeit, ‚in einer Speicherzelle‘‘“ 

Be: arithmetische Operationen mit einem einzigen Befehl ausführen zu lassen. Kombinationen a 

= mit dem technischen Elementarbefehl (5) erlauben, das Ergebnis einer Operation in dem passen- e 

= r den Register für die nächstfolgende bereitzustellen. a 

= | Einige Kombinationen seien an dem. Beispiel der Berechnung von = _ ke) - |d| er- 

läutert. a, b, c und d seien in den Zellen 327 bis 330 gespeichert. Folgende vier Befehle sind‘ a a 
328 p:e N 

E 329 g-+e 
330 bxd 


Das Ergebnis ist nach Zelle 330 geschafft und gleichzeitig zum Abdruck gebracht worden. 


5.4.Besonderheiten bei gleitendem Komma | 


Das Befehlssystem für eine Maschine mit gleitendem Komma läßt sich durch Erweiterung 
deseben angegebenen Grundgerüstes darstellen. Die Speziesoperationen, Verschiebeoperationen 
und u. U. die Sprungbefehle, die im zweiten Code-Verein zusammengefaßt sind, erfahren eine 
Differenzierung. 

Sprungbefehle können sich nunmehr entweder auf die Mantisse oder auf den Exponenten 
beziehen. Diese Alternative kann wiederum durch einen eigenen Digit ausgedrückt werden. 
Die Verschiebeoperationen treten nunmehr auch in Kombination mit einem Herauf- oder 
Herunterzählen des Exponenten auf. Man hat die Varianten 

n-malige Rechtsverschiebung in AC„*) 

q-malige Rechtsverschiebung in AC„, [40,] +g— 40, 


[40,] +g— 46, allein. 


sowie entsprechende für Linksverschiebung. Beschränkt man sich auf nur einmalige Ver- 
schiebungen, so besteht die Möglichkeit, die Verschiebeoperationen mit Speziesoperationen, 
und zwar im unmittelbaren Anschluß an diese, zu kombinieren. 

Dies ist besonders vorteilhaft bei der genauigkeitsrichtigen Ausführung der Multiplikation 
und Division '?), wo abschließend unter Umständen zusätzliche Rechtsverschiebungen mit 
Ausgleichung des Exponenten vorkommen. 


4) ACm bedeutet: Mantissenteil von AC, AC, bedeutet: Exponententeil von AC. 


12) K. Samelson und F. L. Bauer, Optimale Rechengenauigkeit bei gleitendem Komma, 
ZAMP, 4 (1953), S. 312. 


Ai ee eanın En eben Ser ne ee auchi in der von ie 

bei der schwedischen Maschine BARK eingeführten Weise, die sich einfach aus jener bei 
Unterlassung der abschließenden Rechtsverschiebung ergibt, und etwa sogar in einer Aus- 
’führungsart, die die Lage der Mantissen fest läßt). Hier, wie übrigens in diesem ganzen 


Abschnitt, sind Einzelheiten der Durchführung so sehr von technischen Gegebenheiten speziel- 
ler Maschinentypen abhängig, daß wir nicht näher darauf eingehen wollen. Es kam uns nur 
darauf an, die grundlegenden Gedanken in genügender Allgemeinheit zu entwickeln. 


5.5 Einbeziehung der Befehle fürÄnderungsprozesse 
Der Befehl zur Ausführung der Änderung dritter Art ist ein eigener Elementarbefehl, der 
in Kombination mit fast allen bisher aufgeführten Befehlen vorkommen kann. Es mag zweck- 
mäßig sein, für ihn einen eigenen Code-Verein zu gründen, der dann auch die übrigen für 
Änderungen notwendigen Befehle aufnehmen kann. Ein solcher ist der Befehl zur Einstellung 


des Adressen-Addenden-Registers. Dazu gesellt sich der Befehl zur Einleitung der Änderungen 


zweiter Art. Werden zur Kennzeichnung der änderungsbedürftigen Befehle die Vorzeichen- 
stellen von Mantisse und Exponent verwendet, so wird das Befehlssystem nicht weiter belastet. 
Änderungen erster Art sind schließlich schon bei der Eingabe durchzuführen. Sie werden 
also zweckmäßig durch externe Befehle (Betriebszeichen) vorbereitet und ausgelöst. 
Darüber hinaus ist @ie Verwendung solcher Betriebszeichen zum Aufruf einer Befehlssteuerung 
oder eines Teilprogramms für die Konvertierung bereits während der Eingabe ein letzter 
Schritt zur Entlastung des Programms, der sich besonders bei einer dual arbeitenden Maschine 
mit gleitendem Komma empfiehlt. 
Anhang: 
Durchführung der Automatisierung der Änderungsprozesse 
Im folgenden sollen Möglichkeiten der Realisierung der in $3 entwickelten Prin- 
zipien gezeigt werden. Hierfür ist nun ein bestimmter Maschinentyp zugrunde zu legen. Er 
soll sich an das Münchener Projekt anlehnen, das heißt: halblogarithmische Parallelmaschine 
mit einem Befehl pro Wort und Einadreßcode, Gesamtplan etwa wie das folgende Block- 


diagramm !*). 


Bild1 


13) Wenn alle Rechengrößen mit dem Exponenten Null eingegeben werden, bewirkt diese Ausführungsart, 
daß die Maschine wie mit festem Komma rechnet, solange kein Überlaufen eintritt. Kommt dies aber vor, so muß 
die Rechnung noch nicht abgebrochen werden, sondern die Überlaufkorrektur des gleitenden Kommas tritt in 
Kraft. Für Fälle, in denen dies absolut unzulässig ist, wie etwa beim Rechnen mit doppelter Genauigkeit, kann 
durch eine Schaltstellung statt dessen beim Überlauf "Alarm ausgelöst werden. 

14) Herrn W. Proebster danken wir für eingehende Diskussionen über Einzelheiten der technischen 


Durcbführung. 


win neı 

Der normale Eingabevorgang verläuft etwa in: 
(zum Beispiel ein Lochstreifensender) liest die . 
Programmstreifen ab und gibt sie nach dem Akkumu 
gebaut, so wird es auf ein Schlußzeichen nach derjenigen 
jeweiligen Stande des Befehlszählers BZ entspricht. Danach 
weitergeschaltet und die Ablesung des nächsten Wortes beginnt. 
zeichen ‚„‚Anfang eines Teilprogramms‘‘ wird der jeweilige Inhalt des. 
Benützt man hierfür nach dem Vorbild der EDSAC eine feste S 

_ Rechenwerk frei für beliebige Operationen, zum Beispiel für die Umree] 
ins Dualsystem. Das externe Zeichen „Änderung erster Art‘‘ (Bei der EDS ler Kon 
buchstabe ©) bewirkt dann das Aufsuchen dieser Speicherzelle und die Addition ihres Inh 


auf das Befehlswort in AC, das anschließend in die für es selbst bestimmte Speicherzelle trans-- 


feriert wird. Will man aus irgendeinem Grund auf die Adressenumrechnung verzichten ®), 
so kann man sich das zweimalige Ansteuern der festen Speicherzelle sparen und die Bezugs- 
‘adresse aus BZ nach dem Multiplikanden-Register MD schaffen. Die Bestandteile jedes 
Befehlswortes können dann direkt nach AC eingespeist werden, ohne daß MD überhaupt ver- 
wendet wird. Das Zeichen ‚‚Änderung erster Art‘ hat dann nur die Addition von MD auf AC 
zu bewirken. 

B. Änderungen zweiter Art 

Die schrittweise Ausführung eines Programms geht normalerweise nach folgendem 
Schema vor sich: 2 

1. Der Inhalt von BZ geht nach dem Adressen-Register AR. Schalter 1 steht auf I. 

2. AR sucht die Zelle [AR] auf, ihr Inhalt [[AR]] wird in das Befehlsregister BR trans- 
feriert. 

3. Enthält der eben abgelesene Befehl eine Datumadresse, so geht diese, das heißt der In- 
halt des Adressenteils BRA von BR, nach AR. Der Schalter $8 1 wird auf II gestellt. 

4. Die durch den Operationsteil des Befehls BROP gekennzeichnete Operation wird ausge- 
führt. Ist der Befehl ein Sprungbefehl, so wird der Schalter $ 2 auf II gestellt und [BRA] nach 
BZ transferiert. Sonst wird [BZ] um eins erhöht. Anschließend beginnt der Vorgang bei 1. 
von neuem. ; 

Für die Änderungen zweiter Art ist nun dieser Ablauf so zu modifizieren, daß eine feste 
Operation, die Adressensubstitution, an einer Serie von aufeinanderfolgenden Befehlsworten 
entweder ausgeführt wird oder nicht. Als Kennzeichen für die zu ändernden Befehle steht zum 
Beispiel das Mantissenvorzeichen zur Verfügung, als Schlußzeichen des Teilprogramms das 
Exponentenvorzeichen. Bei dem wohl einfachsten Verfahren werden diese beiden Stellen in 
MD abgetastet. Wir wollen dabei das Vorhandensein eines Adressen-Addenden-Registers 
AAR voraussetzen, das als einfacher Zähler arbeiten kann und eine Additionsvorrichtung nach 
AR besitzt. Unbedingt notwendig ist es für den jetzt zu beschreibenden Prozeß zwar nicht, es 
läßt sich jedoch vorteilhaft in ihn einbauen, obwohl es an sich für einen anderen Zweck einge- 
führt wird. 

Der Ablauf der Substitutionsroutine beginnt mit folgendem Zustand: Der aufrufende 
Sprungbefehl ist im Befehlsregister BR eingetroffen. Der Inhalt [BZ] = n des Befehlszählers 
BZ ist die Adresse der Speicherzelle, aus der der Sprungbefehl abgelesen wurde. Der In- 
halt[ BRA]— m des Adressenteils von BR ist die Adresse des ersten Befehls des aufgerufenen Teil- 
programms. Das Adressen-Addenden-Register AAR ist gelöscht. Schalter 8 1 wird auf II ge- 
stellt, Schalter S 2auf I. 

Die Substitutionsroutine läuft dann wie folgt ab: 

1. [BZ] wird um eins au an +»+1(=0,1,2,-- -)erhöht und nach ARgebracht. AR 
sucht die Zelle [AR] auf. Ihr Inhalt, die nächste zu substituierende Adresse, wird nach AC 
transferiert. 

2. [BRA] wird nach AR gebracht und durch Addition von [AAR] auf m + BiR=0T, 
2, ++ )erhöht. Der Inhalt der entsprechenden Zelle, der nächste zu prüfende Befehl des Teil- 
programms, wird nach MD gebracht. 

3. Die Vorzeichenstellen der Mantisse und des Exponenten in MD werden abgetastet: 

Die Kombination 0; 0 bedeutet „keine Substitution, Fortsetzung der Prüfung“. Sie be- 
wirkt Löschung von MD, Erhöhung von [AAR]auf u + 1, Rückkehr nach 2. 

Die Kombination ZL; 0 bedeutet ‚Substitution, Fortsetzung der Prüfung“. Sie bewirkt 
Substitution des Inhalts von AC’in den gerade geprüften Befehl, dessen Adresse noch in AR 


2 15) Adressen sind dann im Octal- oder besser im Sedezimal-System anzugeben. 


3 j 
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j Fre melson u. Bau er, Maßnahmen zur Erzielung us. -. 
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Steht, seh von MD und AC, Erhöhung von [AAR] umeins, Rückkehr nach 1. Eine Lan 
zweiter Stelle bedeutet in jedem Falle „Schluß der Prüfung‘. Sie bewirkt an Stelle der Er- 
höhung von [AAR] und der Rückkehr nach 1. bzw. 2. Ausführung des Sprunges in das Teil- 
programm, das heißt Transfer von 32 = m über Schalter Anach BZ, worauf der Rechen- 
vorgang beginnen kann. 

Es könnte gelegentlich elnais sein, in der Sequenz n -+ v der auf den Sprungbefehl 
folgenden Speicherzellen nicht die Adressen der Rechengrößen, sondern die Rechengrößen selbst er 
unterzubringen'®). Dann ist in den entsprechenden Befehl des Teilprogramms natürichncht 
der Inhalt der Zelle n + », sondern die Adresse n + » selbst zu substituieren. Grundsätzlich RT 
sollten dann beide Möglichkeiten der Substitution offengehalten werden. Dies läßt sich auch > 
durch volle Ausnützung der vier verschiedenen Vorzeichenkombinationen erreichen: Be 


Die Kombinationen 0; 0 und Z; 0 haben dieselbe Bedeutung wie oben, 0; Z bedeutet : 
„Substitution der Adresse n + »v“ en L; L bedeutet ‚‚Schluß‘“ in Kombination mit einer 
festen der drei vorherigen Operationen. er 

Bei diesem komplizierteren Verfahren ist nun zuerst der vorher unter 2. erscheinende Prüf- = 
schritt auszuführen, der aber eben deshalb in AC unter Benützung der dort vorhandenen Ab- 8 
tasteinrichtung ausgeführt werden kann. Die Kombination 0; 0 hat denselben Effekt wie vor- 
her. L;0 bewirkt zunächst Löschung von AC, Transfer von [BZ] nach AR, Aufsuchen der = 
Zelle [AR] und Transfer ihres Inhalts nach AC. 0; L löst stattdessen den Transfer des Inhalts 
i [BZ] von BZ nach AC aus. In beiden Fällen ist nun noch einmal die Adresse des eben geprüften 5 
Befehls durch Transfer von [BRA]nach AR und Addition des noch nicht erhöhten [AA R] her- > 
zustellen. Anschließend wird der Inhalt von AC substituiert. Die Kombination L; Z bewirkt > 


3 
f 
# 


den Sprung in das Teilprogramm durch den Transfer von [B RA] = m nach BZ wie oben. Dieser 
- Sprung kann aber jetzt, wie erwähnt, nur mit einer bestimmten Substitution in den Schluß- 
befehl des Teilprogramms verknüpft werden. Legt man fest, daß der endgültige Ausgang eines 
Teilprogramms immer vomletzten Befehlauserfolgt ‚wassich stets, gegebenenfalls mit Hilfe eines 
zusätzlichen Sprungbefehls, erreichen läßt, so kann man hier zu einer, vollautomatischen Rück- 
kehr in das ‚„‚Rahmenprogramm‘‘ kommen. Unter dem Rahmenprogramm ist dabei diejenige 
Befehlssequenz zu verstehen, die den Sprung (aus Zelle n) in das Teilprogramm enthält. Zu 
diesem Zweck muß die Schlußkombination L; Lerst die durch 0; ZL gekennzeichnete Operation 
und anschließend den Transfer von [BRA]nach BZ auslösen. Dann steht in dem letzten Befehl 
des Teilprogramms automatisch die Adresse des ersten Befehls, der auf die Serie der soeben in 
das Teilprogramm substituierten Größen folgt, und dieser muß nun den Anschluß an die folgen- 
den Programmteile herstellen. 

Andere Möglichkeiten für die Durchführung des Prozesses sind denkbar, bringen jedoch 
nichts wesentliches Neues. Nur auf einen Punkt wollen wir noch hinweisen: Wie eingangs er- 
wähnt, ist das relativ kostspielige Adressen-Addenden-Register AAR für die geschilderte Rou- Be 
tine nicht notwendig. Es genügt vielmehr, wenn der Adressenteil des Befehlsregisters BRA zu 
einem einfachen Zähler ausgebaut wird, den man aber zweckmäßig herunter- und nicht 
heraufzählen läßt. Das Teilprogramm wird dann vom letzten zum ersten Befehl durchgeprüft 
und korrigiert, und beim Eintreffen des Schlußzeichens, das nun den ersten Befehl des Teilpro- 
gramms kennzeichnet, enthält BRA automatisch die für den Beginn der Rechnung benötigte 
Anfangsadresse. Dieses Verfahren erfordert einen ganz geringfügigen technischen Aufwand, 
hat jedoch nur dann einen Sinn, wenn aufein AAR vollständig verzichtet werden soll. 


C. Änderungen dritter Art 


Wie früher erwähnt, sind Änderungen dritter Art von so vorübergehender Natur, daß sie 
zweckmäßig erst dann vorgenommen werden, wenn der betreffende Befehl zur Ausführung in 
das Befehlssteuerwerk gelangt. Hierfür ist zunächst das Adressen-Addenden-Register AAR auf 
die zu addierende Schrittzahl einzustellen. Nach der Überführung der Datumadresse des Be- 
fehls aus BRAnach AR hat der Operationsteildes Befehls durch einen eigenen Digit die Addition 
von [AAR]auf AR zu bewirken. Danach wird die Zelle [AR] aufgesucht und der Befehl ausge- 
führt. Problematisch ist hierbei die Einstellung von AAR. Eine wesentliche Programmver- 
kürzung läßt sich im Falle einer einfachen Induktionsschleife erreichen, in der alle Adressen 
entweder fest bleiben oder in gleicher Weise heraufgezählt werden. Dann kann der ge- 
samte Abzählvorgang in das AAR verlegt und mit der Adressenkorrektur direkt verkoppelt 
werden, und den eigentlich arithmetischen Bestandteilen des Teilprogramms sind insgesamt noch 
drei Befehle hinzuzufügen, nämlich der Befehl, der AAR erstmalig einstellt, und zwei Sprung- 
befehle, von denen der eine die Schleife schließt und der zweite den Ausgang zu den anschlie- 


16) Den Hinweis auf diese Möglichkeit verdanken wir Herrn Dr. R. Piloty. 


x E rn 8 = = 
=. 


Ihr entspricht das Programm | 


0 1— N nach AAR 
1 [ag] nach MR. 
-2 [d&]- [MR] nach AC 
93 TA0]=> 18] mache 
4 Sprung nach 1, wenn AAR <0, [AAR] " ı nach AAR 
5 a nach BR: N 


Be: Sarishlen a ea h. die bei Alert tun nrereh durch die Adressen der Ay Velen 
BT ponenten ay und by ersetzt werden. Der Zählvorgang in AAR wird jeweils durch die Ausführung 
des Sprungbefehls Nr. 4 ausgelöst. AAR enthält daher für jedes n die Zahln— N, undde 
=, effektiven En der Befehle Nr. 1 und 2 sind infolgedessen aa +n— N =a,bzw. 
9 bx + N — =D 
Au Das ee ist leicht zu modifizieren, sobald mehrere Indizes unabhängig abgezählt x 
Ber. werden müssen, also für Systeme von verschränkten Schleifen. In diesem Fall müßte jeder 
E Index in einem eigenen Register abgezählt werden können, und alle diese Register müßten un- 
Eu abhängige Additionen ihrer Inhalte auf AR gestatten — eine technisch durch ein separates 
z Schnellspeicherwerk vor AAR erfüllbare Forderung. 
= 
Ei 
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Wenn man Modifikationen des in Bentwickelten Verfahrens für Änderungen ‚zweiter Art 
or im Sinne des Aufbauprinzips des $ 4 vorsieht, kann man für einfache Serienoperationen, die nur 
> einen Rechenvorgang an zwei Variablenreihen enthalten, zu extrem kurzen Programmen 
5 | kommen. Die notwendigen Modifikationen sind: 

1. Ersetzung der Substitution durch eine arithmetische Operation. 

2. Ersetzung des Prüfschritts im Rechenwerk durch eine Abzählung in AAR wie in dem > 
oben angegebenen Programm. 

3. Streichung des Sprungtransfers von [BRA]nach BZ am Schluß. Der den Prüfvorgang 
einleitende Sprungbefehl wird damit zu einem Serien-Rechenbefehl, dessen Datumadresse die 
letzte Komponente des einen Vektors von Rechengrößen anruft, während der zweite Vektor 
auf den Befehl selbst folgt. Wählt man als arithmetische Operation die Addition,so wird das 
Programm ‚Addition zweier N-komponentiger Vektoren‘‘auf 2 Befehle, die Einstellung von AAR 
und den Serienrechenbefehl, zusammengedrängt. Man könnte sogar den Einstellbefehl noch ein- 
sparen, wenn man die freien Stellen des Rechenbefehls (etwa die Exponentenstellen) für die 
Aufnahme der Komponentenzahl N vorsieht und eine direkte Verbindung zwischen den Expo- 
nentenspuren des Speichers und AAR schafft. 

In derselben Weise kann man auch Operationen ermöglichen, die einen Vektor und eine 
skalare Größe betreffen. Hierfür ist es nur notwendig, den Zählvorgang in BZ abschaltbar zu 
machen, was sich ohne Schwierigkeiten erreichen läßt. Damit ließe sich (unter Ausnützung 

- aller Möglichkeiten) die Multiplikation eines Vektors mit einem Skalar durch einen, das Skalar- 
produkt zweier Vektoren mit zwei Befehlen ausführen, wobei allerdings einer der Ausgangs- 
vektoren verloren geht. 

Vom Standpunkt des Programmierenden aus betrachtet sind diese letzten Verfahren ohne 
Zweifel optimal. Es soll nicht übersehen werden, daß es sich hierbei um eine sehr eng begrenzte 
Klasse von Operationen handelt, deren Automatisierung wohl für Spezialmaschinen von 
großem Nutzen sein kann, aber in gewissem Widerspruch zu dem Gedanken der Allzweck- 
maschine steht. Andererseits aber ist zu beachten, daß eine konsequente Anwendung des 
skizzierten Prinzips in gewisser Weise die Nutzbarmachung der Vorteile eines Zwei-Adreß-Code 
für eine Maschine mit Ein-Adreß-Code darstellt. 


Eingegangen am 20. Juli 53. 


Z. angew. Math. Meclı. 
Bd. 34 Nr.7 Juli 1954 
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KLEINE MITTEILUNGEN 


Lösung spezieller gewöhnlicher Differential- 
gleichungen durch unendliche Operatoren. 
Anwendung auf das ballistische Problem. 


Bei einfach gebauten gewöhnlichen Differentialglei- 
chungen, die indessen keineswegs elementar lösbar sein 
müssen, kann man oft die rekursiv arbeitende Methode 
der sukkzessiven Approximation durch eine indepen- 
dente Darstellung der Lösung ersetzen. Für letztere 
verwenden wir auch die Bezeichnung ‚Unendlicher 
Operator‘ und geben im folgenden einige Beispiele. 

Der unendliche Operator 


67 z z 
yaltffatfsatf.. de...de 0 
ist, falls f(x) beschränkt und stetig im Integrations- 
intervall ist, konvergent, und es gilt 


E7 
SI(a)dz 
EC EN Se (2). 
Der Beweis ergibt sich sofort, wenn man auf die Diffe- 
rentialgleichung y’ =y f(x) mit den Anfangsbedingun- 


genz==(,y = 1die Methode der sukkzessiven Appro- 
ximation anwendet. 


Wir betrachten jetzt die Differentialgleichung 


in der (x, y)-Ebene, von der ein schmaler Streifen längs 
y =0 ausgeschlossen sei. Die Anfangsbedingungen 
seien =(0,y=% = 0, und zur Abkürzung schreiben 
wir 


z 
F («) er fie) de. 
Die Lösung von (3) ist 
Er 
d 
y-4u+Fr+ ou (4). 
Yt F+ z 
„tfF+ Nah 
0 0 
0 
(Yy +9) 
Die Differentialgleichung 
y’=ya)/y+f@), =0, y=0":(d), 


hat für y > 0 die Lösung 


zT 
vor [vVr+fviwHnaie 
0 


= 
hierbei ist F (a2) = f(x) de. 
0 


In (3) sei f, in (5) f und-y iin dem betrachteten Inter- 
vall beschränkt und stetig vorausgesetzt. 

Um einen Anhaltspunkt für das Konvergenztempo 
der unendlichen Operatoren zu gewinnen, betrachten 


wir beispielsweise (5) und (6) für f(z) = 0. Man gewinnt 
die Identität 


Tz 
1 Maren 
27 yY(z) de = fvY Fornäede (7). 
P 0 0 


Setzt man in (7) y(z) =2, so ergibt die schrittweise 
Auswertung 


E b 
= .lim'a,r" 
no 


mit 
4 
a ee) 
bn 
n=1,2,3. 
Nach fünf Schritten wird a, = 1,092, b, = = ; nach 
: 1023 
zehn Schritten @10= 1,005, 54, = 1094° 


Aus (1) und (2) erhält man für f(x) =1 die Darstel- 
lung 
z 3 zT 
e=ltfir fit fit...de 
0 
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In (8) erreicht man nach n Schritten 


Zum Schluß machen wir noch eine Bemerkung zum 
Hauptproblem der ebenen Ballistik. 

Bedeuten v die Geschoßgeschwindigkeit, f(v) die 
Funktion des Luftwiderstandes, ce den Geschoßfaktor, 
ö die Neigung der Babntangente gegen die Horizontale, 
so gelten die bekannten Differentialgleichungen 


v= — cfwW) —gsind .....(9 
ven EEE, (92): 
vv Rest 


Bild1 


V=V* 


Bild 2 


Durch Elimination von ® bekommt man 
d + ef)? + wN? = g? . (10), 
durch Differentiation von (9°) nach der Zeit ergibt sich 
bei Benützung von (9°) 
dv 


—vı=—cf (v) od + vr 
dv 
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und hieraus bei Beachtung von (10) 
d 
dv 
oder Een f , 
wi = 20 (ed wit 2e@F + 2) 
® 


u. 


Für den gesamten aufsteigenden, sowie für den ab- 
steigenden Ast der Flugbahn bis zum Minimalwert der 
Geschwindigkeit v*, für den verschwindet, ist die Wur- 
zel mit positivem, für den Rest der Bahn mit negativem 
Vorzeichen zu nehmen. Den qualitativen Verlauf von 
(v ö) und von (v ö)? deuten die beiden Bilder an. 


Die Pfeilrichtungen auf den Kurven entsprechen dem 


zeitlichen Ablauf. +. 5 

' Die Differentialgleichung (12) für (v 0)? als ‚Funktion 
von v ist von der gleichen Gestalt wie (5), die Lösung 
also durch einen Operator der Form (6) darstellbar !). 


BUCHBESPRECHUNGEN Be: 


'Tables des Fonctions de Legendre associ6es. 


Fonction.associde de premiere espece P} (cos Ö), premier 
fascieule [n: —0,5 (0,1) 10; m: 0 (1) 5; 9: 0 (1) 90°]. 
Caleuldes pour le Centre National d’Etudes des 
Tel&communications (C.N.E.T.) par le Service 
M&canographique de l’Armöe de Terre. (Collection 
technique et scientifique du C.N.E.T.) XXII + 
292 S. Paris 1952. Editions de la Revue d’Optique. 


Der Anlaß zur Berechnung der vorliegenden Tafel der 
zugeordneten Legendreschen Funktionen erster 
Art D* (cos 9) — der noch weitere folgen sollen — 


war ein Problem aus der Praxis: Bei der numerischen 
Berechnung von Doppelkonusantennen benötigt man 


Werte von 20 (cos ®) für reelles, nicht-ganzzahliges n. 


Bisher existierten jedoch nur Tafeln für ganzzahlige 
Werte von n. 


Der Tafel wurde die Ho bson sche Definition 
d” P„ (cos ©) 

d (cos 9)” 
zugrunde gelegt, die sich von der sonst üblichen (vgl. 
z.B. Jahnke-Emde) durch den Faktor (— 1)” unter- 
scheidet. Für die Bereehnung wurde der m-n-Bereich 
in Felder unterteilt, deren jedes einen Wert von m 
und zehn Werte vonn,n=r-+0,1s(r ganz, s= (0, 
l,... 9), umfaßt. Für zwei Felder -— m=4r=5 
und m =5,r =5 — wurde zunächst 

m 
Dr (cos d) 

sin” 9 
aus der Reihenentwicklung mit großer Genauigkeit 
direkt ermittelt. Die Berechnung von 7% (9) in den 
übrigen Feldern wurde dann mit Hilfe von Rekursions- 
formeln vom allgemeinen Typ 


AFTER LCeTEN=-0 


2 (cos #) = (— 1)” sin® d 


une 


a,b,c,d ganz 

ATBACHE team. 
durchgeführt. Aus den so errechneten Werten von 
7? (8) wurde durch Multiplikation mit sin” ® 


schließlich Pr (cos 9) bestimmt. Die Zahl der in der 


Tafel angegebenen geltenden Ziffern der Funktions- 
werte bewegt sich — bedingt durch die Benutzung der 
Rekursionsformeln — zwischen 3 und 11. Den Zahlen- 
tafeln sind noch Höhenkarten der Funktionen 


PR} (cos 9) (m= 0, 1, 2, 3, 4, 5) beigegeben, die einen 
schnellen Überblick über deren Verlauf gestatten. 


„Le eften "(iD | 
V® w 


en Kurven: 

Verzweigungspun! 
den Zweige bat sein 
ginge aus von den Anfangswerte 
Variierte man dann v* und c, 
alle Schußbahnen. Um die zu 

€, vg, 9, gehörende herauszufinden, hät 
vo, = — ef(w) —gsind, 


zu setzen und dann diejenige Schußbahn auszusuchen, 1 
bei-der im 1. Teil ö, und v, zusammengehören. 
Konstanz. b "R. ‚Müller. 


ee ee Le \ # 
1) In (6) wäre # durch y, + F zu ersetzen und die Integrationen von 
z, bis x zu erstrecken. 


Außer bei dem oben genannten speziellen Problem . = 
wird die Tafel z. B. bei der numerischen Behandlung en - 


von partiellen Differentialgleichungen, bei denen die 
Randbedingungen auf Rotationsflächen 2. Ordnung 
gegeben sind, Verwendung finden. 


Dresden. A. Schubert. 


Dr.-Ing.R.Zurmühl,PraktischeMathema- 
tik für IngenTeure und Physuker 
X +481 Seiten mit 114 Abb. Berlin/Göttingen/Hei- 
delberg 1953. Springer-Verlag. Preis geb. 28,50 DM. 

Das vorliegende Buch gibt eine moderne Darstel- 
lung der für den Ingenieur wichtigsten Methoden der 
praktischen Mathematik. Es dient der Erweiterung 
des für den Hocbschulingenieur in der Praxis wich- 
tigen mathematischen Rüstzeug. Entsprechend 
dieser Zielsetzung werden keine langen Beweise durch- 
geführt. Der Verfasser legt vielmehr Wert darauf, 
den mathematischen Kern der Probleme deutlich zu 
machen und die grundlegenden Ideen ausführlich zu 
erörtern. Das ist ausgezeichnet gelungen. Der größte 
Teil des Buches ist den linearen Problemen gewidmet, 
der Verfasser verwendet dabei durchweg die Matrizen- 
schreibweise, wodurch die Darstellung sehr übersicht- 
lich wird. Jeder Abschnitt enthält mehrere voll- 
ständig durchgerechnete Beispiele, die das Einarbei- 
ten in, die vorgeführten Methoden sehr erleichtern. 
Der Verfasser berücksichtigt die neueste Entwicklung, 
z. B. sind das Summenverfahren nach Sassen- 
feld als Erweiterung des Mehrstellenverfahrens 
und die Newtonsche Verbesserung von Matrizen- 
eigenwerten nach Unger wohl erstmals in ein Lehr- 
buch aufgenommen. 

Im I. Kapitel wird Zahlenrechnen mit Rechen- 
schieber und Rechenmaschine behandelt. Der Re- 
ferent würde es begrüßen, wenn in einer Neuauflage 
in dieses Kapitel einiges über Nomographie aufgenom- 
men werden könnte. Das II. Kapitel bringt Iterations- 
verfahren zur Lösung von Gleichungen mit einer Un- 
bekannten, Verfahren zur Lösung algebraischer Glei- 
chungen und Stabilitätskriterien. Im III. Kapitel 
wird zunächst der Gaußsche Algorithmus in ursprüng- 
licher und mechanisierter Form elementar behandelt. 
Danach wird ein kurzer Abriß des Matrizenkalküls 
gegeben, worauf die Matrizenschreibweise kon- 
sequent verwendet wird. Es folgen die Hessenberg- 
sche Aufstellung der charakteristischen Gleichung 
und die Iterationsverfahren. Wenn auch dem In- 
genieur die Darstellung in Matrizenform oft noch un- 
geläufig ist, so kann er sich doch an Hand der Para- 
graphen über Matrizenrechnung und der zahlreichen 
Beispiele jederzeit die Bedeutung wieder ins Gedächt- 
nis rufen. Im IV. Kapitel folgen Interpolation und 


.- - 


Buchbesprechungen 


+ 


Integration, wobei auch die Lagrange-Hermiteschen 
Interpolationsformeln berücksichtigt sind. V. und 


‚VI. Kapitel bringen ausführlich Ausgleichsrechnung 


und Darstellung willkürlicher Funktionen. Im VII.Ka- 
pitel werden Anfangswertprobleme bei gewöhnlichen 
Differentialgleichungen behandelt, im VIII. Kapitel 
die Rand- und Eigenwertaufgaben vor allem bei line- 
aren Differentialgleichungen. In diesen letzten beiden 


“3 Kapiteln sowie im III. und IV. Kapitel sind die mo- 


dernen Verfahren besonders stark berücksichtigt. Auf 
die Behandlung von partiellen Differentialgleichungen 


“und Integralgleichungen wird verzichtet, da dieses den 


Rahmen des Buches gesprengt hätte und ja hierfür 
auch genügend Spezialliteratur existiert. 

Das Werk wird vor allem auch dem schon länger in 
der Praxis stehenden Ingenieur und Physiker von 
Nutzen sein können, da es überall an die modernsten 
Methoden der praktischen Mathematik in leicht ver- 
ständlicher Weise heranführt. Es ist ihm eine recht 
große Verbreitung in Kreisen der Physiker und In- 
genieure zu wünschen. Der Verlag hat das Buch in 
bekannter Weise erstklassig ausgestattet. 


Dresden. Bachmann. 


Dr. J. E. Hofmann (Honorarprof. a. d. Univ. Tü- 


bingen, Geschichte der Mathematik 


Erster Teil: Von den Anfängen bis zum Auftreten 
von Fermat und Descartes (Sammlung Göschen 
Bd. 226). 200 S. Berlin 1953. Verlag W. de Gruyter 
& Co. Preis 2,40 DM. 

Die vorliegende außerordentlich gedrängte Dar- 
stellung der Geschichte der Mathematik bis ein- 
schließlich der Zeit des Frühbarock geht auch auf die 
Geschichte der Nachbargebiete, insbesondere auf die 
geistigen Strömungen des Mittelalters ein. Die antike 
Mathematik wird ganz kurz behandelt. Das Bänd- 
chen bringt eine erdrückende Fülle von Namen und 
Material, doch geben ausführliche Quellennachweise 
am Ende jeden Abschnittes und ein Namen- und 
Schriftenverzeichnis von 43 Seiten Umfang am Schluß 
des Bandes Wegweiser für das eingehendere Studium 
von Einzelheiten. 

Dresden. Willers. 

Heinz Prüfer, Projektive Geometrie. 
Aus dem Nachlaß herausgegeben von G. Fleddermann 
und @. Köthe. 2. unveränd. Aufl. VII + 314 S. m. 
3 ganzseitigen Abb. und 251 Fig. im Text. Leipzig 
1935. Akademische Verlagsgesellschaft Geest u. 
Portig K. G. Preis-geb. 9,— DM. 

Nach dem Tode von Prüfer (1934) gaben G. Fledder- 
mann und G. Köthe 1939 sein Manuskript mit den 
notwendigen Ergänzungen heraus. Davon liegt jetzt 
eine 2., im wesentlichen unveränderte Auflage vor. 

Die projektive Geometrie hat sich heute gegen zwei 
übermächtige Nachbarn zu verteidigen, die sie zu er- 
drücken drohen: die Axiomatik und Grundlagen- 
forschung und die lineare Algebra. Prüfer versteht es, 
beiden gegenüber das Eigenrecht der projektiven 
Geometrie zu wahren, sie beide als Hilfsmittel zu ge- 
brauchen, ohne sie überwuchern zu lassen. Er ver- 
steht es in sehr schöner Weise, auf der einen Seite die 
großen Gesichtspunkte herauszuarbeiten, auf der 
anderen uns eine Fülle von reizvollen konkreten geo- 
metrischen Einzelheiten vorzuführen. 

Er beginnt mit einer ausführlichen Begründung 
durch die Axiome der Verknüpfung, der Anordnung 
und der Stetigkeit. Gleich zu Anfang führt er die un- 
eigentlichen Elemente und die beiden für die projek- 
tive Geometrie wesentlichen Übertragungsprinzipe 
ein: das Dualitätsprinzip und das der Invarianz 
gegenüber Projektionen. Mit ihrer Hilfe begründet 
er die Lehre von den projektiven Abbildungen syn- 
thetisch. Die metrische Geometrie wird durch neue 
Axiome in die projektive Geometrie eingebaut. Das 
Verhältnis kehrt sich aber dann um, wenn er den 


Fundamentalsatz der projektiven Geometrie nun aus 
den Kongruenzaxiomen, ohne Benutzung der Stetig- 
keit ableitet. 

In ganz natürlicher. Weise schließt sich die Einfüh- 
rung der nicht-euklidischen Geometrie an. Verf. geht 


dabei ganz konkret von Konstruktionen in einem be- 


grenzten Raumstück aus und gewinnt unter Benut- 
zung der Axiome der Orthogonalität und Stetigkeit 
die drei Formen der nicht-euklidischen Geometrie. 
Von dem gewonnenen Standpunkt aus erschließen 
sich die’ Prinzipien der darstellenden Geometrie in 
sehr einleuchtender Weise. Verf. führt sie durch bis 


‘zur Behandlung von Schattenkonstruktionen, Durch- 


dringung von Tetraedern, Tangenten an Regelfläche 
und Reliefperspektive. i 

Erst in letzten Kapitel führt er Koordinaten ein. 

Das Buch ist sehr sorgfältig bis in alle Einzelheiten 
durchdacht. Es führt den Anfänger in geschickter 
Weise in die Tatsachen und Denkweise des Gegen- 
standes ein. Es kann sowohl in Verbindung mit einer 
Vorlesung als auch im Selbststudium dank seiner 
klaren und leicht verständlichen Darstellung mit Vor- 
teil benutzt werden. Aber je weiter man fortgeschrit- 
ten ist, um so mehr bemerkt man erst den eigentlichen 
Tiefgang des Werkes; es vermittelt die Einsicht in die 
Zusammenhänge und in die kristallene Ordnung der 
Tatsachen und führt die Beweise stets in Einordnung 
unter große Gesichtspunkte und unter Hinweis auf 
das Wesentliche. Auch der reifste Mathematiker wird 
aus diesem Buch noch lernen können und an ihm seine 
Freude haben. 


Halle. Ott-Heinrich Keller. 


.B. H. Jennings (Professor am Technologischen In- 
stitut der Northwestern University) und W. L. Rogers 
(Assistent Professor am Technologischen Institut der 
Northwestern University) Gas Turbine Ana- 
lysis. and Praktice. VII -+487S. m, zahl- 
reichen Figuren. New York/Toronto/London 1953. 
Me Graw-Hill Book Company. Preis geb. 64 sh. 

Die Autoren haben ihr Buch in der Absicht ge- 
schrieben, eine ausführliche Darstellung über die 
Arbeitsweise und das Anwendungsgebiet der Gas- 
turbinen zu geben. Dabei wird die Gasturbine im 
Zusammenhang mit anderen Wärmekraftmaschinen 
gesehen und die gemeinsamen Grundlagen in der 
Mechanik und Thermodynamik werden eingehend 
klargelegt. 

Zu begrüßen sind die für die Berechnung des Ar- 

eitsprozesses notwendigen Enthalpie—Entropie-Ta- 
feln von Luft, denen spektroskopische Messungen der 
spezifischen Wärme zugrunde gelegt sind. 

Im Kapitel über die Entwurfsgrundlagen werden 
u.a. Ratschläge zum Entwurf der Schaufelprofile 
in der konventionellen Art und Weise erteilt, wie es 
beim Bau von Turbinen und Pumpen üblich ist. Da- 
bei werden die Strömungsvorgänge nur mittels der 
sogenannten Stromfadentheorie berechnet. Jedoch 
sollte die Gitterströmung zweidimensional behandelt 
werden. Erst damit erhält der gestaltende Ingenieur 
den notwendigen Einblick in die Einzelheiten des 
Strömungsvorganges und lernt die Güte eines Pro- 
files zu beurteilen. 

Erfreulich ist wiederum das Eingehen auf zahlreiche 
Probleme. Neben konstruktiven Hinweisen, z.B. 
bei Schwingungsfragen sowie bei der Bebandlung der 
Festigkeitsrechnung unter Berücksichtigung der Wär- 
mespannungen, werden die Werkstofffragen mit vielen 
Einzelheiten über die Wärmefestigkeit besprochen. 
Das Kapitel über die Brennstoffe belehrt an Hand von 
Tabellenmaterial über physikalische und chemische 
Eigenschaften wichtiger Treibmittel. 

Das Buch beschäftigt sich in den Anwendungen 
vorwiegend mit stationär eingebauten Gasturbinen, 
jedoch wird auch auf andere Bauarten, z. B. bei Loko- 
motiven, eingegangen und den Flugzeugtriebwerken 


19* 


Propellerantrieb mit Abgasturbine gegenüber ge- 
stellt wird. ‚ } 


Wenn die Verfasser die Beratung des praktischen . 


Ingenieurs und die Einführung von Studierenden in 
dieses Spezialgebiet beabsichtigten, wird das Buch 
seinen Zweck erfüllen. 


Dresden. W. Albring. 


Uno Olsson, Non-cireular eylindrical 
gears. XI + 216 S. Stockholm 1953. Acta Poly- 
technica. Mechanical Engineering Series Vol. 2 
Nr. 10. Preis 18,— Schw. Kr. 

Es gibt Maschinen und Geräte, bei denen die trei- 
benden Wellen gleichförmig umlaufen, während sich 
die getriebenen Wellen mit ungleichförmiger Winkel- 
geschwindigkeit drehen. Diese verlangte Ungleich- 
förmigkeit der Abtriebswelle tritt besonders bei Werk- 
zeugmaschinen, Textilmaschinen und feinmecha- 
nischen Geräten auf und kann entweder durch Ge- 
lenkmechanismen (z. B. Antiparallelkurbel) oder 
durch unrunde (in besonderen Fällen elliptische) 
Zahnräder erreicht werden. Wenn keine Qualitäts- 
anforderungen vorliegen, kann man z.B. elliptische 
Zahnräder in einem einfachen Teilverfahren mit 
Zahnformfräsern anfertigen. Man muß dann mit 
Hilfe der Legendre schen elliptischen Funktionen 
und elliptischer Integrale für jede zu fräsende Zahn- 
lücke die Ordinaten und Einstellwinkel berechnen. 

Das vorliegende Werk will helfen, die Herstellung 
unrunder Zahnräder in kontinuierlichem Fräsver- 
fahren zu ermöglichen und eine hohe Genauigkeit zu 
erreichen. 

Nach einer Einführung im Kapitel 1 werden die 
hauptsächlichsten Gesetze und Konstruktionen der 
Rollkurven für Zahnräder behandelt. Es folgen im 
Kapitel 2 einige Definitionen und die Bedingungen, 
‚die von zwei Kurven gefordert werden, wenn sie ein 
Rollkurvenpaar bilden. 

Im Kapitel 3 werden Gleichungen und Bogen- 
längen von Kurven zweiten Grades behandelt: a) die 
Ellipse, b) die Hyperbel, c) die Parabel, d) die Ellipsen 
mit imaginären Achsen, e) die imaginäre Hyperbel, 
f) die Grenzkurve im besonderen Fall, wenn die Ex- 
zentrizität e = sin a gegen 0 geht. 

Kapitel 4 behandelt die Rollkurven zu den Kegel- 
schnitt-Kurven, insbesondere der Ellipse, Hyperbel, 
Parabel, der imaginären Ellipse und der imaginären 
Hyperbel. 

Kapitel 5 befaßt sich mit der Konstruktion von 
Rollkurvenpaaren für unrunde Zahnräder mittels 
Kegelschnitt-Rollkurven und mit dem Ersatz einer 
gegebenen Rollkurve durch eine Serie von Bögen 
von Kegelschnitt-Rollkurven. 

Kapitel 6 bringt die Herstellung der Werkzeuge, 
insbesondere der Räderfräser für unrunde Zahnräder. 
Um ein reibungsloses gegenseitiges Abwälzen der 
Zahnkurven des Werkzeuges und des Rades zu ge- 
währleisten, ist es notwendig, daß die Geschwindig- 
keit der relativen Bewegung im Berührungspunkt 
der Zahnflanken gleich Null ist. 

Im Schluß-Kapitel 7 wird gezeigt, daß auch Kegel- 
radgetriebe mit unrunden Zahnrädern möglich sind. 

Den Anhang des Buches bilden zahlreiche Tabellen 
und Diagramme (Seite 168— 214), die dem Konstruk- 
teur bei der Bestimmung der Getriebekenngrößen 
einen erheblichen Teil der Rechenarbeit ersparen. 
Mit einem ausführlichen Literaturverzeichnis schließt 
das Buch ab. 

Es ist zu hoffen, daß durch das vorliegende Werk 
die Herstellgenauigkeit der unrunden Zahnräder 
wesentlich gesteigert werden kann. Die Möglichkeit 
der Anwendung solcher Getriebe wird den Industrie- 
zweigen des Maschinen- und Gerätebaues sehr will- 
kommen sein. 


Dresden. W. Liehtenmhertd:, 


ein besonderes Kapitel gewidmet, worin dann die. u 
reine Gasturbine u.a. dem Raketenantrieb und dem 


Nach allger 
der Grundbegr 
lichen Methoden der Err 
und Beschleunigung in eber 
Form entwickelt und auf „Hebel-Gelen 
Kurventriebe und Rädertriebe angewanı - 
sonders ausführlicher Darstellung der letzteren ein- 
schließlich der Entwicklung der Zahnprofile undihrer 
Korrektur. Diesem 262 Seiten umfassenden kine- 
matischen Hauptteil folgt der kürzere Teil über die 
Dynamik der Getriebe. Dieser beschäftigt sich mit 
Reibungskräften, mit den statischen Kräften in Ge- 


trieben und in e großen Zügen mit der Dyna- . 
mik des zentrischen Schubkurbeltriebes, dem Massen- “ 
ausgleich bei drehenden und hin- und hergehenden 
Massen sowie der Regelung der Maschinen. Es wird Er 
also hierbei nicht eigentlich die Dynamik der Getriebe, Es 
sondern nur die Dynamik des zentrischen Schubkurbel- sr 


getriebes geboten und zwar in der üblichen, im deut- - 
schen Schrifttum behandelten Form. Fe 


Das Buch ist bei erwähnter Breite klar und anschau- 
lich geschrieben, es kann als Einführung indie v 
Kinematik zum Selbststudium dienen. Die Schrift- 3 
tumshinweise beziehen sich ausschließlich auf russi- 
sche Arbeiten, abgesehen von wenigen Ausnahmen, 
doch ist es wertvoll, hierdurch einen Einblick in das 
russische Schrifttum zu bekommen. Daß der Druck- 
fehlerteufel Newton in das 20te Jahrhundert 
versetzt, wird der Leser mit Vergnügen bemerken. 


Aachen. W. Meyer zur Capellen. 


Dr. Bernhard Sticker (Prof. a. d. Univ. Bonn), 
FintsitielligemıTafele der u Trigonpe 
metrischen Funktionen. Ausgabe A für 
Altgradteilung. 52 S. + 1 Interpolationstafel. Bonn 
1954. Ferd. Dümmlers Verlag. Preis geb. 7,60 DM. 


Die vorliegende fünfstellige Tafel enthält die nu- 
merischen Werte der sechs trigonometrischen Funk- 
tionen sin, cosec, tg, ctg, sec und cos für jede Bogen- 
minute des Quadranten. Bis auf den Bereich 
0° <as 2° (für cosece und ctg) bzw. 88° <a < 90° 
(für sec und tg) ist durchgängig lineare Interpola- 
tion erlaubt. Zu diesem Zweck sind die ersten Diffe- 
renzen angegeben, falls sie mehr als 5 Einheiten der 
letzten Stelle betragen. Die Interpolation in den 
oben bezeichneten Bereichen wird durch die Angabe 
der Hilfsgrößen 0 =’ -coseca und T=a’.ctga 
— anstelle der Differenzen — sehr erleichtert. 


Auf einen kleinen Schönheitsfehler sei noch hin- 
gewiesen: Im Vorwort wird gesagt, alle Funktions- 
werte seien mit fünf Dezimalstellen oder mit wenig- 
stens fünf geltenden Ziffern angegeben. Dieses Prin- 
zipist an einer Stelle durchbrochen: zwischen 35’ und 
2° sind für cosec«& und ctg& nur vier geltende - 
Ziffern zu finden. 


Dresden. A.Schubert. 


Dr. S. Valentiner (Prof. emer. a.d. Bergakademie 
Clausthal. Vektoranalysis. (Sammlung 
Göschen Bd. 354.) Neudruck der 7. Aufl. 138 S. m. 
19 Abb. Berlin 1954. Verlag Walter de Gruyter. 
Preis geb. 2,40 DM. 


Der vorliegende Neudruck der 7. Aufl. dieses Bänd- 
chens weist gegenüber früheren Auflagen, z. B. der 5., 
die. in Bd. 18 (1938) S. 381 dieser Zeitschrift be- 
sprochen wurde, nur unwesentliche Änderungen und 
Ergänzungen auf. Wie schon dort betont wurde, 
stellt das Bändchen eine gute Einführung in das 
Gebiet der Vektoranalysis dar. 


Dresden. Willers. 
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R. Beyer, Kine matische Getriebe- 


 synthese. (Grundlagen einer quantitativen Ge- 


triebelehre ebener Getriebe.) VIIT+2178. mit 
258 Abb. Berlin/Göttingen/Heidelberg 1953. Sprin- 
ger-Verlag. Preis geb. 36,— DM. 

Wie der Haupttitel des Buches, man möchte sagen 
doppelt ausdrückt, liegt hier die Betonung auf der 
rein kinematischen, bewegungsgeometrischen Seite 
der Getriebelehre, aber immer in Blick auf die tech- 
nischen Anwendungen, d..h. auf den Entwurf von Ge- 
trieben zur Ableitung bestimmter Bewegungen. Um 
unter möglichst einheitlichen Gesichtspunkten dar- 
stellen zu können, werden u.a. infinitesimal benach- 
barte Lagen durch Grenzübergang aus endlich be- 
nachbarten Lagen gewonnen, und so ergibt sich die 
Gliederung: Zwei Lagen eines Getriebegliedes, drei 
endlich, drei infinitesimal benachbarte Lagen, vier 
endlich benachbarte Lagen, Lagenzuordnungen, vier 
und fünf infinitesimal benachbarte Lagen, das Gelenk- 
viereck als Kurbelgetriebe, Ellipsen als Schmiegungs- 
kurven, vorgeschriebene Geschwindigkeits- und Be- 
schleunigungsverhältnisse bei Kurbeltrieben und 
ebenen Kurvenscheiben. 

Das Buch ist klar geschrieben bei ausgezeichneter 
Unterstützung durch gute Abbildungen und eilt bei 
rein theoretischen Erörterungen schnell wieder den 
praktischen Anwendungen zu. Der Verfasser stützt 
sich auf zahlreiche, insbesondere auch neuere Arbeiten, 
wie ein ausführliches, ziemlich vollständiges Schrift- 
tumsverzeichnis am Schluß des Buches und zahlreiche 
Schrifttumshinweise im Text zeigen, stützt sich aber 
auch auf eine Reihe eigener Arbeiten. Erfreulich ist 
hierbei die Beachtung der Arbeiten vonReinhold 
Müller. 

Das ‚für den Konstrukteur, für die Vorlesung und 
das Selbststudium‘ geschriebene und sehr gut aus- 
gestattete Buch bedeutet eine wertvolle Bereicherung 
des getriebetechnischen Schrifttums (es wird auch den 
Mathematiker interessieren), bedarf aber für den Kon- 
strukteur wie auch für den Studenten beim Selbst- 
studium bereits eines gewissen Vertrautseins mit 
kinematischen Dingen. Dem Buch ist weite Ver- 
breitung zu wünschen. 


Aachen. W. Meyer zur Capellen. 


R.Kraus, Grundlagen des systema- 
tischen Getriebeaufbaus. 3798. m. 
243 Abb. Berlin 1952. Verlag Technik. Preis geb. 
19,— DM. 

Der durch seine zahlreichen Veröffentlichungen aus 
dem Gebiet der Getriebetechnik bekannte Verfasser 
— eigenwillig, wie so manche Wissenschaftler in die- 
sem Fach — geht auch in diesem Buch seine eigenen 
Wege, sei es in den Methoden, sei es in der Bezeich- 
nungsweise.Wenn man sich jedoch an die Darstellungs- 
weise gewöhnt hat, d. h. sich ‚‚eingelesen“ hat, ist man 
überrascht, wieviel Anregungen und Hinweise auf 
vielfältige, für den Konstrukteur wichtige Getriebe- 
formen es enthält — selbst wenn man in Einzelfällen 
zur Darstellung und Entwicklung andere Wege gehen 
möchte. 

Das Buch ist in 39 Abschnitte unterteilt. Es be- 
ginnt mit den Begriffen Freiheitsgrad, Beweglich- 
keitsgrad und der Zahlsynthese ebener und räum- 
licher Getriebe. Anschließend werden aus bestimmten 
Mechanismen (die hier noch nicht als Getriebe an- 
gesehen werden) eine Reihe einfacher und zusammen- 
gesetzter Getriebe entwickelt, bei denen naturgemäß 
die Gelenkgetriebe im Vordergrund stehen, doch die 
Räder- und Kurventriebe nicht darüber vernach- 
ässigt werden. Zum Schluß wird noch kurz auf Ge- 
triebe mit Überschluß und auf die Sperr- sowie Rei- 
bungsgetriebe eingegangen. Im Vordergrund steht, 
dem Titel entsprechend, der Aufbau der Getriebe in 
systematischer Weise unter Annahme bestimmter 
Aufgaben des Einzelgetriebes. So treten notwendiger- 


h re Er . 
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. handelten Methoden zurück und wird die*tGetriebe- 
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weise die rein bewegungsgeometrischen, in anderen 
Schriften (auch denen des Verfassers) ausführlich be- 


analyse nur soweit behandelt, wie es zum Verständnis 
notwendig ist. Daher sind auch keine Schrifttums- 
hinweise gegeben und nur neun Namen erwähnt. 
Das Buch wird dem Getriebefachmann willkommen 
on und dem Konstrukteur manche Anregungen ge- 
en. z nr 


Aachen. W. Meyer ZAURT, Capellen. 


“_Dietrieh Küchemann und Johanna Weber, Aero- = 
dynamics of Propulsion. IX.u. 340 S. mit T 
169 Abb. u. 17 Tab. London, New York, Toronto 1953. 


McGraw-Hill Publishing Company Ltd. Preis 64 sh6d. 
Die hydro- und aerodynamische Literatur der ee 
ersten Jahrhunderthälfte stand wesentlich unter dem ee 
Zeichen der Grenzschichttheorie und, soweit sie das Br 
Fliegen betraf, unter dem der Tragflügeltheorie. Beide 2 Er: 
Gebiete haben vielfach zusammenfassende Darstel- Br. 
lungen in Lehrbüchern und anderweitig gefunden, wo- | 


bei der Schwerpunkt der Tragflügeltheorie in der Er- E 
fassung des Auftriebs lag. Nun steht aber die”tech- 3 
nische Entwicklung der letzten zehn Jahre wesentlich a 
unter dem Zeichen neuer Formen des Vortriebs, und f 
unter diesem hat sich ein neues Gebiet auch der Theorie j 
gebildet, die ‚Aerodynamik des Vortriebs“. Wenn sie . 
auch wegen der Vielgestaltigkeit der neuen Formen K 
der Triebwerke keine ebenso einheitliche Darstellung 
gestattet wie die erstgenannten Gebiete, so verlangt - 
doch die technische Entwicklung dringend nach 
einem zusammenfassenden Werk, das der neueren 
Entwicklung Rechnung trägt und auch .die neue 
Literatur entsprechend verwertet. Wissenschaftler 
und Ingenieure werden sehr erfreut sein, in dem vor- 
liegenden Werk diesem Bedürfnis weitgehend Rech- 
nung getragen zu finden. 

Soweit es sich um die grundlegenden theoretischen 
Fragen aerodynamischer und thermodynamischer 
Natur handelt, die alle Formen betreffen, sind sie in 
verschiedenen Kapiteln geschickt zusammengefaßt. 
Darüber hinaus gelten besondere Kapitel jeweils den j 
einzelnen besonderen Triebwerksarten: der” Mantel- ie 
schraube, dem Staustrahlrohr, der Strahlturbine und 
den Wirkungen der verschiedenen Einbauarten der 
Strahltriebwerke. Ein weiteres Kapitel beschäftigt sich 
sich mit Fragen des Strahles und seiner Wirkung auf 
die Nachbarschaft. Es folgt ein biologischer Exkurs 
über das Zusammenwirken in der Erzeugung von Auf- 
trieb und Vortrieb bei fliegenden Tieren. Den Be- 
schluß bildet ein Kapitel über Kühlungsfragen mit ge- 
schiekten Hinweisen für den Konstrukteur und ein 
Anhang mit Funktionstafeln. 

Da das Buch wirklich eine Lücke ausfüllt, wäre es 
sehr zu begrüßen, wenn es bald ins Deutsche zurück- 
übersetzt oder in der wohl vorhandenen deutschen 
Urschrift veröffentlicht würde. 


Weilburg. 


a 


L. Schiller. 


Dr. Fritz Stüssi (o. Prof. a. d. ETH Zürich) Vor- 
lesungen überBaustatik. Bd. 1: Statisch 
bestimmte Systeme — Spannungsberechnung — 
Elastische Formänderungen — Stabilitätsprobleme 
— Seile. 2. Aufl. 370 S. m. 336 Abb. Basel und Stutt- 
gart 1953. Verlag Birkhäuser. Preis geb. 40,05 S Fr., 
brosch. 35,90 S Fr. 

Die\vorliegende 2. Auflage des Buches Baustatik I 
von Stüssi bedarf eigentlich keiner besonderen Be- 
sprechung, dafür ist der Verfasser und das Buch zu 
bekannt und eingeführt. Trotzdem kann nicht ein- 
dringlich genug auf die Vorzüge des Buches hinge- 
wiesen werden. 

Unser heutiges technisches Schrifttum krankt 
daran, daß die einzelnen Veröffentlichungen kaum 
noch lesbar in des Wortes ursprünglicher Be- 


ja 


‚zifferte Aneinanderreihung von Tatsachen dar, bei 


. denen der rote Faden — nicht nur für den Anfänger 


so wichtig zur Einführung in die Problematik und 
‚Methodik der betreffenden Disziplin, sondern ebenso 
für den erfahrenen Praktiker erwünscht —nur schwer 
zu erkennen ist oder überhaupt fehlt. 

Das Buch von Stüssi ist einerühmliche Ausnahme. 
Von der ersten bis zur letzten Seite in einer klaren 
und überzeugenden Form geschrieben, behandelt es 
die statiseh bestimmten Systeme, die klassische 
Biegungslehre, die elastischen Formänderungen, die 
Stabilitätsprobleme und die Statik der Seile. Trotz 
dieses umfassenden Stoffes, der in den Abschnitten 
über Biegung und Stabilität auch die neuesten Er- 
kenntnisse einbeschließt, ist es kein dickleibiger 
Wälzer geworden, sondern ein wirkliches Lehrbuch 
geblieben, das von jedem Studenten bewältigt werden 
kann. 


Braunschweig. Henn. 


Dr. E. Schmidt (o. Prof. a. d. Techn. Hochschule 
München). Einführung in die Techni- 
sche Thermodynamik und in die 
Grundlagen der Chemischen Ther- 
modynamik. 5.ber. Aufl. XVI + 5208. mit 
224 Abb. und 69 Tabellen sowie 3 Dampftafeln 
als Anlage. Berlin/Göttingen/Heidelberg 1953. Sprin- 
ger-Verlag. Preis geb. 30),— DM. 

Die Neuauflage ist bis auf geringfügige Ande- 
rungen, die im wesentlichen die neuen Vereinbarungen 
in der Festlegung der internationalen Temperatur- 
skala auf der 9. Generalkonferenz für Maß und Ge- 
wicht in Paris 1948 betreffen, ein unveränderter 
Nachdruck der vorhergehenden. Es gelten daher für 
die Neuauflage die gleichen Anmerkungen wie für 
jene (s. ausführliches Referat in Z. angew. Math. 
Mech. 31 (1951) Nr. 1/2.) Neben der vorzüg- 
lichen Ausgestaltung der beiden Hauptsätze und 
ihrer technischen Anwendungen auf Gase und 


Dämpfe, der Behandlung der Verbrennung und der. 


Strömungsvorgänge mit ihrer Anwendung auf die 
Strömungsmaschinen, Raketen und Strahlantriebe 
erscheinen die Kapitel über die Wärmeübertragung 
und die Strahlungsübertragung, in denen z. B. die so 
wichtige Gasstrahlung nur gestreift ist, und besonders 
das kaum 2 Seiten umfassende Kapitel über die 
Diffusion allzu dürftig. Dafür wären die seit der 
4. Auflage hinzugekommenen Kapitel XXIund XXII 
entbehrlich gewesen, die sich mit Teilgebieten der 
chemischen Thermodynamik, nämlich der Wärme- 
tönung, dem Massenwirkungsgesetz und dem ther- 
mischen Wärmesatz beschäftigen. Für den Ma- 
schinen- und Wärmeingenieur sind sie bedeutungslos, 
belasten zwecklos das sonst so wertvolle Werk, das 
zweifellos zu den besten deutschen Lehrbüchern der 
technischen Thermodynamik zählt. 


Dresden. H. Ra] tim 


Dr. E. Trost, (Prof. a. Technikum Winterthur) 
Primzahlen (Rlemente der Mathematik vom 
höheren Standpunkt aus. Bd. II). 95 S. Basel und 
Stuttgart 1953. Verlag Birkhäuser. Preis brosch. 
13,50 SFr. 

Auf 95 Seiten gibt der Verf. eine nach systema- 
tischen Gesichtspunkten geordnete Übersicht über 
die wesentlichen Ergebnisse der Primzahltheorie. 
Trotz des beschränkten Umfanges werden zahlen- 
theoretische Kenntnisse nicht vorausgesetzt, sondern 
in den beiden ersten Abschnitten in dem erforder- 
lichen Umfange entwickelt. 

Die Darstellung ist von gedrängter Kürze. Den- 
noch dient das in der Sammlung ‚Elemente der 
Mathematik vom höheren Standpunkt aus“ er- 
schienene Bändchen nicht bloß der Orientierung 
durch Nachschlagen. Drei Umstände machen die 


der wichtigsten Ve; J 


‚deutung sind. Sie bieten sich meist nur als eine be- I 


ebenfalls nur zu begrüßen. 


Dresden. Draeger. { h 


Review of Input and Output Hanıpt 
mentusedin Computing Systems (Joint 
AIEE-IRE-ACM Computer Conference) 142 S. m. 
zahlr. Abb. New York 1953. American Institute 
of Electrical Engineers. Preisgeh. 4$. 5 

In der Zeit vom 10.-bis 12. Dezember 1952 fand in 
New York eine von rund 1100 Personen besuchte Ta- 
gung statt,in der technische Fragen der Ein- und Aus- 
gabevorrichtungen von Ziffernrechenautomaten be- 
handelt wurden. In vorliegendem Heft sind die Vor- 
träge und Diskussionsbeiträge dieser Tagung wieder- 
gegeben. Zunächst wurden die für die Ein- und Aus- 
gabe zur Verfügungstehenden technischen Möglichkei- 
ten allgemein besprochen. Dabei handelt es sich in 
der Hauptsache um Einrichtungen für den Übergang 
von geringer zu hoher Geschwindigkeit bei der Ein- 
gabe und umgekehrt bei der Ausgabe. Weiterhin 
wurden dann im einzelnen diese Einrichtungen beim 
SEAC, UNIVAC, RAYDAC und IBM 701 erörtert. 
Der Rest der Vorträge beschäftigt sich im wesent- 
lichen mit den verschiedenen Ausgabevorrichtungen. 
Aber auch andere Themen werden behandelt wie der 
Übergang vom Analogie- zum Ziffernrechner, nu- 
merisch gesteuerte Mahlvorrichtungen usw. Ins- 
gesamt erhält man einen Überblick über den Stand 
der Entwicklung dieses Gebietes Ende 1952. 


Dresden. Willers. 


R.v. Mises (Prof. a. d. Harvard Universität Cam- 
bridge) u. Th.v. Kärmän (Prof. a. d. Columbia Uni- 
versität, New York), Advances in Applied 
Mechanics. Band II. %X+324S. m. zahlr. 
Abb. New York 1953. Academic Press Inc. Preis 
geb. 9,— $. 

Wie in den bisher erschienenen Bänden dieser Buch- 
reihe wird in einzelnen Aufsätzen über den gegen- 
wärtigen Stand der Forschung auf dem Gebiete der 
angewandten Mechanik in Amerika und England 
berichtet. Im ersten Abschnitt gibt G F. Carrier 
einen Überblick über spezielle Randwertaufgaben, 
wobei Probleme der Ozeanströmung, der Wärme- 
übertragung, der Konvektion, der Relaxations- 
schwingungen des Van der Pol'schen Oszillators u.a. 
behandelt werden. Noch bestehende Schwierigkeiten 
und Nachteile der einzelnen Integrationsverfahren 
werden erläutert. ©. Zaldastani behandelt die 
eindimensionale isentropische Flüssigkeitsströmung. 
Durch Einführung neuer Variablen ergibt sich eine 
übersichtliche Darstellung der allgemeinen Lösung, 
welche für verschiedene Parameterwerte, insbesondere 
Wellenbewegungen, sowie Strömungsvorgänge im ge- 
schlossenen Rohr näher untersucht wird. F.N. 
Frenkiel behandelt ein Problem der turbulenten 
Diffussion, und zwar die Verteilung der Haupt- 
konzentration in einem Strömungsfeld mit homo- 
gener Turbulenz. Mit den üblichen Definitionen und 
entsprechend den Methoden der Statistik werden die 
Grundgleichungen für turbulente Diffusion aufgestellt 
und für den Fall einer Punktquelle bei verschiedenen 
Bedingungen näher diskutiert. H.F. Ludloff 
behandelt Stoßvorgänge in der Aerodynamik. Die 
Auslösung von Schallwellen wird näher untersucht 
und führt auf Erweiterungen des Hu yghensschen 
Prinzips. G. Guderley betrachtet Stoßwellen bei 


gemischten Unterschall-Überschallströmungen, wo- 
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bei Annahmen zugrunde gelegt werden, welche eine 
verhältnismäßig übersichtliche Darstellung der End- 
ergebnisse ermöglichen. L. Rosenhead behan- 
delt Wirbelsysteme im Strömungsschatten von Kör- 
pern und zeigt, daß die Folgerungen aus den bisherigen: 
Formulierungen der Wirbelströmung an gewisse Ge- 
nauigkeitsgrenzen gebunden sind. H. Geiringer 


berichtet über neuere Ergebnisse zur Theorie des 


idealplastischen Körpers. Nach allgemeinen Dar- 
stellungen des dreidimensionalen und des ebenen 
Problems wird die Lösung für das ebene Problem mit 


Hilfe des Charakteristiken-Verfahrens im Sinne be- 


kannter Darstellungen von Mises, Neuber und 
Sauer angegeben. Nach Formulierung des voll- 
ständigen Problems (einschlie ßlich der Verschiebungen) 
werden Wellenbewegungen untersucht und Beispiele 
diskutiert. A.I. Bellin betrachtet nichtautonome 
Systeme, d.h. nichtlineare Differentialgleichungen 
zweiter Ordnung, bei welchen die zweite Ableitung 
der gesuchten Funktion x nach der Zeit nicht nur 
von & selbst und seiner ersten Ableitung, sondern 
noch explizit von der Zeit abhängt. Mit Hilfe einer 
topologischen Transformation läßt sich eine Klassifi- 
zierung der möglichen Systeme durchführen. Mit den 
einzelnen Bedingungen lassen sich Aussagen über die 
Stabilität des jeweiligen Lösungstyps verbinden. 


Dresden. H. Neuber. 


Alwin Walther (Prof. a. d. T. H. Darmstadt). 
Angewandte Mathematik Tell u II. 
(Naturforschung und Medizin in Deutschland 1939 
bis 1946, Bd.3 u. 4. Für Deutschland bestimmte 
Ausgabe der Fiat Review of German Science ) Wein- 
heim/Bergstraße 1953, Verlag Chemie GmbH. Teil I: 
401 S. m. 113 Abb. Preis kart. 20,— DM; Teıl I: 
165 S. m. 22 Abb. Preis kart. 10,— DM. 

Von den beiden als letzte erschienenen Bänden der 
Fiat-Berichte ist der erste numerischen, instrumen- 
tellen und graphischen Methoden der angewandten 
Mathematik gewidmet. In zwölf verschieden langen 
Abschnitten werden behandelt: Graphische und 
numerische Verfahren (Collatz), Praxis der kon- 


formen Abbildung (Ullrich), Nomographie und 
Rechenschieber (Walther und Kron), Mathe- 
matische Maschinen und Instrumente. Instrumentelle 
Verfahren (Walther und Dreyer), Mathe- 
matische Zahlentafeln. Numerische Untersuchung 
spezieller Funktionen (Walther und Unger), 
Wahrscheinlichkeitsrechnung und Mathematische 
Statistik (Schulz), Versicherungs-, Wirtschafts- 
und Finanzmathematik (Lore y), Anwendung der 
Mathematik auf Biologie, Medizin und Bevölkerungs- 


- wissenschaft (Geppert) Mathematische und ma- 


schinelle Methoden beim Chiffrieren und Dechiffrieren 
(Rohrbach), Darstellende und konstruktive 
Geometrie (Rehbock), Vektor- und Tensor- 
rechnung (Lyra), Verfolgungskurven (Hose- 


.mann). Im zweiten Band finden sich die wichtigsten 


Teile der Mechanik und zwar: Theoretische Mechanik 
(Hamel), Allgemeine Elastizitätstheorie (Schmidt), 
Baustatik (Bornscheuer), Statik im Schiffbau 
(Schnadel), Spannungsoptik (Föppl und 
Mönch), Schwingungen mit endlich vielen Frei- 
heitsgraden (Klotter), Schwingungen elastischer 
Körper (Schmidt), Flügelflattern (Küssner), 
Flugmechanik (Stüper), Theorie der Regelung und 
Steuerung (O p pelt), Schwingungsmessung (Kl ot- 
ter), Fahrmechanik (Maier und Riekert) 
Plastizitätstheorie (Pöschl). Die meisten Referate 
enthalten ausführliche Schrifttumsangaben. Ihre 
Abfassung hat zum Teil mehr oder weniger große 
Schwierigkeiten gemacht, da die erforderlichen Unter- 
lagen vielfach vernichtet waren. Trotzdem ist es ge- 
lungen, einen umfassenden Überblick der wissen- 
schaftlichen Entwicklung der betreffenden Gebiete 
in Deutschland in den Jahren 1939—1946 zu geben, 
der, anschließend an den Stand der vorhergehenden 
Zeit, zusammen mit der nun wieder zugänglichen Ent- 
wicklung im Ausland als Grundlage für weitere Ar- 
beiten dienen wird. Herausgeber wie Verfassern der 
einzelnen Aufsätze muß man für die wohl oft unter 
nicht leichten Verhältnissen geleistete Arbeit herzlichen 
Dank aussprechen. 


Dresden. Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen, ausführliche Besprechung bleibt vorbehalten. 


Prof. Dr. A. Lottermoser (ehem. Direktor des 
Kolloidehemischen Instituts an d. Techn. Hochschule 
Dresden). Kurze Einführung in die 
Kolloidcehemie unter besonderer 
Berücksichtigung der anorgani- 
schen Kolloide. Dritte, erg. Auflage von 
Dr.-Ing. Carl Kalauch. (Forschungsinst. f. chem. 
Technologie Meinsberg bei Waldheim). VIII + 260 S. 
m. 111 Abb. Dresden/Leipzig 1954. Verlag von 
Theodor Steinkopff. Preis geb. 10,— DM. 


Keith Swainger, Ph. D. (Research Fellow in the 
University of London). Analysis of Defor- 
mation. ‚VokT. Mathematical Theory. 
XIX +285 S. London 1954. Chapman & Hall Ltd. 
Preis geb. 63. 


G.Pickertt, Lineare Algebra — Nor- 
malformen von Matrizen. (Enzyklo- 
pädie der Mathematischen Wissenschaften mit Ein- 
schluß ihrer Anwendungen, Band I 1, Heft 3, Teil I, 
6/7.) 728. Leipzig 1953. B. G. Teubner Verlags- 
gesellschaft. Preis geb. 7,50 DM. 


Kurt Reidemeister (o. Prof. a. d. Universität Mar- 
burg). Topologie der Polyeder und 
Kombinatorische Topologie der 
Komplexe. (Mathematik und ihre Anwendungen 
in Physik und Technik. Reihe A. Band 17). Zweite 


Aufl. XI + 196 S. m. 69 Abb. Leipzig 1953. Akade- 
mische Verlagsgesellschaft Geest & Portig K.G. 
Preis geb. 14,— DM. 


W.Maak, Darstellungstheorie un- 
endlicher Gruppen und fastperio- 
dische Funktionen. (Enzyklopädie der 
Mathematischen Wissenschaften mit Einschluß ihrer 
Anwendungen, Band Il, Heft 7, Teil I, 16.) 26 S. 
Leipzig 1953. B.G. Teubner Verlagsgesellschaft. 
Preis geb. 3,— DM. 


Paul F. Byrd and Morris D. Friedman, Hand- 
bookofelliptiece Integrals for Engi- 
neers and Physicists. (Die Grundlehren 
der Mathematischen Wissenschaften in Einzeldar- 
stellungen Bd. LXVII.) XIII + 355 S. mit 22 Abb. 
Berlin / Göttingen / Heidelberg 1954. Springer-Verlag. 
Preis geh. 36,— DM, geb. 39,60 DM. 


Dozent Ove Petterson, Circular Plates 
Subjectedto Radially Symmetrical 
Transverse Load Combined with 
Uniform Compression or Tensionin 
the Plane of the Plate. (Tekniska Skrifter 
utgivna av Svenska Teknologföreningen nr 153) 
30 S. m. 21 Abb. Stockholm 1953. Teknisk Tidskrifts 
Förlag. Preis geh. 7 kr. 
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Dr.-Ing. Rudolf Richter (Prof.em.a.d. T. H. Karls- 
ruhe), Elektrische Maschinen. Bd. Ill. 
Die Transformatoren. 2. Aufl. 321 8. m. 
230 Abb. Basel/Stuttgart 1954. Verlag Birkhäuser. 
Preis geb. 26schw. Fr. (26,— DM). 


Prof.Dr. F.H. Müller, Das Relaxations- 
verhalten der Materie. (2. Marburger 
Diskussionstagung) IV + 224 S.m. 122 Abb. Darm- 
stadt 1953. Verlag Dr. Dietrich Steinkopff. Preis 
brosch. 24,— DM. 


Dr.-Ing.B.Eck, Technische Strömungs- 
lehre. 4. verb. Aufl. X +4228S. m. 407 Abb. 
Berlin/ Göttingen/Heidelberg 1954. Springer-Verlag. 
Preis geb. 29,40 DM. 


Nachrichten — Zuschriften an den Herausgeber 


Z. angew. Math. Mech. 
* Bd.34 Nr.7 Juli 1954 


Dr.-Ing. Kurt Beyer (ehem. Prof. a. d. T. H. Dres- 
den, Technische Mechanik für Bau- 
ingenieure. (Herausgegeben von Dipl.-Ing. 
H. Baldauf und Dipl.-Ing. C.-H. Hänig.) 
VIII +1598S. m. 326 Abb. Leipzig 1954. Verlag 
S. Hirzel. Preis geb. 11,80 DM. 


Dr.-Ing. M. Schuler (Prof. em. der Universität Göt- 
tingen,, Einführung in die Mechanik. 
Teill. Mechanik des Massenpunktes. 
202 S. m. 114 Abb. Teil I. Mechanik-der 
Punktsysteme. 2108. m. 166 Abb. Wolfen- 
büttel 1950/51. Verlag Benno Kracke. Preis geb. 
je 12,— DM. 


NACHRICHTEN ° 


Berlin. Am 15. Januar verstarb auf der Fahrt 
von Berlin nach Dresden unerwartet der Professor 
an der Technischen Universität Berlin-Charlotten- 
burg Dr.-Ing. HermannAlt. 


München: In Gemeinschaft mit der Firma Zeiß- 
Aerotopograph München veranstaltete das Institut 
für Photogrammetrie, Topographie und Allgemeine 
Kartographie vom 15.— 27. 3. 54 einen internationalen 
Kurs in Photogrammetrie, der unter dem Namen 

„Münchner Photogrammetrische 
Wochen 1954“ die Tradition dieser früher in 
Jena abgehaltenen Lehrgänge fortsetzte. Die Kurs- 


leitung lag in den Händen der Herren Professoren 
Dr. R. Finsterwalder und Dr.K. Schwi- 
defsky. 


München. Mit Wirkung vom 15. Februar 1954 
wurde bei der Technischen Hochschule München eine 
Dokumentationsstelle auf dem Gebiet der Luftfahrt- 
forschung errichtet. Sie führt die Bezeichnung ‚‚Doku- 
mentationsstelle auf dem Gebiet der Luftfahrtfor- 
schung, insbesondere auf dem Gebiet der Elektro- 
technik, der Flugmeteorologie, des Flugzeugbaues 
und der Triebwerkforschung bei der Technischen 
Hochschule München‘, 


ZUSCHRIFIEN AN DEN HERAUSGEBER 


Zu F.W.Schäfke, Verbesserte Konver- 
genz- und Fehlerabschätzungen für die Stö- 
rungsrechnung, Z. angew. Math. Mech. 33 
(1953), S. 255— 259, 

Am Schluß seiner interessanten Abhandlung ver- 
gleicht Herr F. W. Schäfke verschiedene untere 
Schranken 0, bis 0, für den Konvergenzradius der 


Potenzreihen (*) & Amr uk sowie die zugehörigen Ab- 
k—=1 


l 
schätzungen des-Fehlers (£) Am (1) — Am — & Ammpuk 


(z. B. 0,<< 0,, siehe letzte Formelzeile, dort befindet 
sich übrigens ein Druckfehler, es soll offenbar 
80, = 2% = 0, statt 80, = 2%= 0, heißen). Die 
dortige Darstellung könnte vielleicht zu Mißverständ- 
nissen Anlaß geben. Deshalb sei dazu hier folgendes 
bemerkt. 

Mit o, bezeichnet Herr Schäfke wohl die in meiner 
Dissertation genannte Schranke 0,,, = [2a +b +e]7, 
worin (mit den Bezeichnungen von Herrn Schäfke) 


1 
a=y - lach Eml) und-c= = (Gym, Ym)| ist und 
dm dm Inn 


für b entweder y (2 u oder = \@yml| gesetzt 
m dm dm 

werden kann. Diese in den von mir behandelten 
Fällen allgemein verwendbare Schranke kann nicht 
größer als o,ausfallen. Fürhermiteschen Stör- 
operator jedoch wurde in meiner Dissertation (S. 13) 
eine verbesserte Schranke 09,5 —= [a +2b +c]! an- 
gegeben. Diese kann größer als 0, werden, wie das fol- 
gende (als Grenzfall zu wertende) Beispiel zeigt: 


F Are 7)y=39, m=3, 


hier ist Gym = 0, bei Verwendung einer Funktion 
= const-& erhält man daher 9,» = 20,1 = 2 


Die Abschätzungen des Fehlers (#) wurden in 
meiner Dissertation unter Benutzung von Majoranten 
für die Reiben (*) hergeleitet. Bei hermiteschem Stör- 
operator ist der Konvergenzradius der verwendeten 
Majorante (mindestens) 0,,,, bei nichthermiteschem G 
jedoch nicht o,,,, sondern (mindestens) 


093 = [la +2fadb+c]!. 
Auch o,,; kann größer als 0, werden. Beispiel: 


— y’tucos2nz  y’=4Ay, yV=y(l)=0, 


n“ 
Im= — 
Mm 4 ’ 
man erhält bei Verwendung von y= n und b= Em 


die Zahlen 0,3 = a so — Es gibt daher auch 


Fälle — sowohl solche mit hermiteschem als auch 
solche mit nichthermiteschem Störoperator —, in 
denen die Abschätzungen von Herrn Schäfke für 
den Fehler ($) bei „gro ße m“ I nicht genauer ausfal- 
len als die mit Hilfe der erwähnten Majoranten 
hergeleiteten. 


Zu obigem möchte ich jedoch abschließend bemer- 
ken, daß in der Zusammenfassung der wesentlichsten 


‘ Ergebnisse am Schluß meiner Dissertation die Grö- 


Ben 03, und E,,, nicht erwähnt sind, da es mir (neben 
der Schematisierung der Ergebnisse) hauptsächlich 
auf die Verbesserung der Abschätzungen des Fehlers 
(*) für kleines } ankam. 


Hannover. J. Schroeder. 
Die vorstehenden Ausführungen haben Herrn 
Dr.Schäfke vorgelegen. 
Willers. 
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